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Model  laws  governing  the  design  of  qeumechar.ical  model  studies  of  underground  open¬ 
ings  in  rock  subjected  to  static  loads  are  developed  using  dimensional  analysis  and  the 
theory' of  models.  The  significant  variables  influencing  the  behavior  of  rock  considered 
in  this  study  are  the  free-field  stresses,  the  properties  of  the  discontinuities  in  the 
rock  mass,  the  opening  geometry,  and  the  structural  liner  supporting  the  opening.  The 
prototype  chosen  for  study  is  a  short  section  of  a  long  circular  tunnel  which  is  under¬ 
ground  at  a  depth  of  more  than  four  tunnel  diameters. 

The  development  of  geomechanical  modeling  techniques  for  the  construction  and  test¬ 
ing  of  jointed  rock  models  Is  described.  The  strength  properties  of  the  intact  model 
material  and  of  the  Joint  surfaces  are  also  presented.  The  models  were  tested  in  plane 
strain  (zero  strain  parallel  to  the  tur.nel  axis)  by  the  use  of  a  controlled  rigid  longi¬ 
tudinal  loading  head. 

Five  jointed  models  were  tested  to  determine  the  effect  of  tunnel  liner  stiffness 
and  the  ratio  of  joint  spacing  to  tunnel  diameter  on  the  behavior  of  turlnel  liners  in 
jointed  rock  masses.  The  model  tunnels  were  lined  with  plexiglass  liners  of  varying 
thickness  to  simulate  reinforced  concrete  liners  In  a  jointed  rock  mass.  The  effect  of 
the  ratio  of  the  tunnel  diameter  to  joirt  spacing  on  structural  behavior  was  studied  by 
testing  three  different  size  tunnels  in  the  model  reck  mass  composed  of  1  - 1 n .  square 
joint  blocks.  The  structural  behavior  of  the  mode’s  was  assessed  by  means  of  quanti¬ 
tative  measurements  such  as  measurement  of  the  diameter  changes  of  the  liner,  the 
extensometer  measurements  of  displacements  in  the  rock  mass  behind  the  tunnel  wall. 

The  observed  model  behavior  was  analyzed  and  compared  to  a  previously  developed 
elasto-plastic  analysis.  The  analyses  showed  that  elastic  theory  was  sufficiently 
accurate  for  calculating  the  diametrical  strains  of  the  lined  tunnels  up  to  diametrical 
strains  of  about  1.2X.  For  diametrical  strains  above  1.2X,  the  elasto-plastic  theory 
was  used  successfully  to  calculate  the  tunnel  diametrical  strains  within  an  error  of 
about  ♦  15X. 
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ABSTRACT 


Model  laws  governing  the  design  of  gcon.  'i  hanical  model  studies  of  underground 
openings  in  rock  subjected  to  static  loads  are  developed  using  dimensional  analysis  and 
the  theory  of  models.  The  significant  variables  influencing  the  behavior  of  rock  con¬ 
sidered  in  this  study  arc  the  frcc-ficld  stresses,  the  properties  of  the  discontinuities  in  the 
rock  mass,  the  opening  geometry,  and  the  structural  li  ter  supporting  the  opening.  The 
prototype  chosen  for  studv  is  a  short  section  of  a  lo’  g  ciicular  tunnel  which  is  under¬ 
ground  at  a  depth  of  more  than  four  tunnel  diamete  s. 

The  development  of  geomechanical  modeling  techniques  for  tire  construction  and 
testing  of  jointed  rock  models  is  described.  The  strength  properties  of  the  int-'Ct  model 
material  and  of  the  joint  surfaces  are  also  presented.  The  models  were  tested  in  plane 
strain  (zero  strain  parallel  to  the  tunnel  axis)  by  the  use  of  a  controlled  rigid  longitudinal 
loading  head. 

I’ ive  jointed  models  were  tested  to  determine  the  effect  of  tunnel  liner  stiffness  and 
the  ratio  of  joint  spacing  to  tunnel  diameter  on  the  behavior  of  tunnel  liners  in  jointed 
rock  masses.  The  model  tunnels  were  lined  with  plexiglass  liners  of  varying  thickness  to 
simulate  reinforced  concrete  liners  in  a  jointed  rock  mass.  The  effect  of  the  ratio  of  the 
tunnel  diameter  to  joint  spacing  on  structural  behavior  was  studied  by  testing  three 
different  size  tunnels  in  the  model  rock  mass  composed  of  1-in.  square  joint  blocks.  The 
structural  behavior  of  the  models  was  assessed  by  means  of  quantitative  measurements 
such  as  measurement  of  the  diameter  changes  of  the  liner,  the  extensometcr  measure¬ 
ments  of  displacements  in  the  rock  mass  behind  the  tunnel  wall. 

The  observed  model  behavioi  was  analyzed  and  compared  to  a  previously  developed 
clastoplastic  analysis.  The  analyses  showed  that  clastic  theory  was  sufficiently  accurate 
for  calculating  the  diametrical  strains  of  the  lined  tunnels  up  to  diametrical  strains  of 
about  1.2%.  For  diametrical  strains  above  1.2%,  the  clastoplastic  theory  was  used  suc¬ 
cessfully  to  calculate  the  tunnel  diametrical  strains  within  an  error  of  about  ±  15%. 
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MODEL  TESTS  OF  LINED  TUNNELS 
IN  A  JOINTED  ROCK  MASS 


1  INTRODUCTION 

Justification  for  Geomechanical  Model  Studies.  The 

ament  volume  nt  underground  construction  activity 
tor  military  and  civil  work  is  unprecedented  in  history, 
and  is  expected  to  expand  significantly  in  the  future. 
Moreover,  the  size  and  complexity  of  underground 
structures  is  steadily  increasing.  Although  significant 
advances  have  been  made  in  recent  years,  our  under¬ 
standing  of  the  behavior  of  underground  openings  is 
pom  because  the  ptoblem  is  extremely  complex  and 
highly  indeterminate  due  to  the  discontinuous  nature 
of  jointed  rock  masses.  A  more  comprehensive  under¬ 
standing  of  the  interaction  of  various  structural  liners 
with  jointed  rock  masses  is  necessary  if  future  under¬ 
ground  works  are  to  be  designed  and  constructed 
rationally  and  economically. 

Among  the  variables  which  exert  a  significant  in¬ 
fluence  on  the  behavior  of  underground  openings  in 
rock  are  the  following:  natural  free- field  stresses  in 
rock  mass;  artificial,  superimposed  loads  coming  cither 
from  within  the  opening  or  applied  to  the  rock  mass  at 
some  distant  point;  mechanical  properties  of  the  intact 
rock  material;  the  nature  of  discontinuities  in  the  rock 
mass,  such  as  joints,  faults,  and  bedding  planes;  the 
geometry  of  the  opening;  the  structural  liner  or  rock 
bolts  supporting  the  opening;  and  the  techniques  and 
sequence  hv  which  the  opening  is  constructed. 

Analytical  methods  of  predicting  the  behavior  of 
underground  openings  in  rock  are  quite  limited  in  ap¬ 
plicability  because  the  jointed  rock  mass  is  generally 
discontinuous;  the  construction  procedures  arc  of  great 
significance  but  arc  often  indeterminate;  and  the  sys¬ 
tem  consisting  of  the  underground  structure  and  the 
surrounding  rock  mass  is  also  highly  indeterminate. 
Solutions  from  the  theory  of  elasticity  arc  directly 
applicable  only  to  a  limited  number  of  rock  masses 
whose  properties  approach  the  assumptions  of  clastic 
theory  because  joint  spacings  are  large  and  the  stress 
levels  imposed  arc  below  failure  stresses  for  the  rock 
mass. 

Numerical  analyses  utilizing  finite  clement  tech¬ 
niques  and  electronic  computers  are  much  more  versa¬ 
tile  and  are  becoming  quite  sophisticated.  They  arc 


capable  of  considering  to  varying  degrees  such  proper¬ 
ties  as  openings  of  any  shape,  anisotropy,  non- homog¬ 
eneous  layering,  non-linear  clastic-plastic  frictional 
material  properties,  and  discontinuities  in  the  rock 
mass  (see  for  example  Reyes.1  Goodman,2  and  Zien- 
kicwicz.3).  At  present,  however,  restrictions  of  com¬ 
puter  size,  development  of  finite  element  techniques, 
and  our  ability  to  accurately  determine  and  describe 
the  actual  complex  properties  of  the  intact  rock  ma¬ 
terial,  (lie  rock  mass,  and  the  underground  structures 
limit  the  usefulness  of  this  method  of  analysis. 

Hence,  in  addition  to  the  use  of  theoretical  and 
analytical  techniques  in  studying  and  predicting  the 
behavior  of  underground  openings,  it  is  necessary  to 
use  empirical  techniques  in  field  and  laboratory  experi¬ 
mentation.  Limited  field  data  concerning  the  behavior 
of  underground  openings  subjected  to  static  and  dy¬ 
namic  loadings  arc  available,  but  it  is  not  possible  to 
extrapolate  these  data  directly  to  predict  the  behavior 
at  other  sites  where  the  pertinent  variables  have  differ¬ 
ent  values.  In  order  to  do  this  one  must  have  a  quanti¬ 
tative  basis  for  determining  how  changes  in  the  vari¬ 
ables  will  influence  the  behavior  of  the  opening.  To 
develop  empirically  such  quantitative  relationships 
between  the  pertinent  variables  and  the  behavior  of  the 
opening,  data  must  be  obtained  ovr,  i  wide  range  of 
the  variables.  The  cost  and  impractical  ty  of  obtaining 
data  from  many  full-scale  field  construction  sites  limits 
the  usefulness  of  this  approach.  Field  data  does,  how¬ 
ever.  offer  the  only  ready  means  of  studying  the  influ¬ 
ence  of  construction  techniques.  In  addition,  concepts 
and  predictions  derived  from  other  methods  of  study 
must  ultimately  be  checked  by  field  observations  to 
substantiate  or  disprove  their  validity  and  usefulness. 


'  S.E.  Reyes.  Elastic  Plastic  Analysis  of  Underground  Open¬ 
ings  by  the  Finite  Element  Method,  Ph.  D.  Thesis  (Universi¬ 
ty  of  Illinois,  1 966). 

>  R.K.  Goodman.  “On  the  Distribution  of  Stresses  Around 
Circular  Tunnels  in  Non-Homoguious  Rocks."  Proceedings. 
1st  International  Congress,  Vo!  2  (lnt.  Soe.  of  Rock  Mec.. 
1966)  pp  249-255:  R.  E.  Goodman.  R.L.  Taylor  and  T.  L. 
Hrckke.  “A  Model  for  the  Mechanics  of  Jointed  Rock." 
Proceedings  ASCII.  Vol  94.  No.  I  (May  1968)  pp  637-659. 
’  O.O.  Zienkiewicz.  “Continuum  Mechanics  as  an  Approach  to 
Rock  Mass  Problems, "  Rock  Mechanics  in  Engineering  Prac¬ 
tice,  Stagp  and  Zienkiewicz,  cd.  (John  Wiley  &  Sons.  1968). 
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A  most  promising  technique  lor  studying  the  in- 
11  nonce  ot  many  oi  the  pertinent  variables  appears  to 
he  the  use  ol  geomechanieal  moilels.  In  this  technique, 
a  small  scale  tnorlel  ol  the  underground  opening  is  con¬ 
structed  m  a  material  which  accurately  models  the 
propel  lies  ol  the  actual  rock  mass  in  the  field.  The 
model  is  then  loaded  m  such  a  manner  as  to  reproduce 
the  stress  state  which  exists  in  the  real  prototype 
underground.  If  the  requirements  of  similitude  are 
satisfied,  the  behavior  of  the  model  then  reproduces 
the  behavior  of  the  prototype  in  all  respects:  distri¬ 
bution  of  stresses,  distribution  of  strains  and  defor¬ 
mations  both  elastic  and  inelastic,  and  failure  modes. 
There  are  relatively  few  theoretical  limitations  on  this 
technique.  Its  usefulness  appears  to  depend  mostly 
upon  the  skill  and  ingenuity  of  the  investigator  in  solv¬ 
ing  the  many  practical  and  technical  problems  involved 
in  trying  to  satisfy  the  requirements  of  similitude,  par- 
ticularly  in  modeling  the  details  of  the  underground 
structure  and  of  the  geologic  environment  of  the  proto¬ 
type.  These  problems  are  not  insignificant. 

The  validity  and  tisefulne.s  of  structural  models  is 
well  established  in  many  phases  o.  engineering  research 
and  design.  The  most  notable  examples  in  civil  engi¬ 
neering  probably  are  the  structural  model  tests  of  arch 
dams  and  other  structures  at  such  places  as  Laboratorio 
Nacional  de  Engcnharia  Civil  in  Lisbon.  Portugal  and 
Istituto  Spcrimentalc  Modelli  c  Strutturc  in  Bergamo, 
Italy.  The  next  step  beyond  structural  modeling  is  the 
use  of  geomechanical  models  in  which  not  only  the 
proposed  engineering  structure  is  modeled,  but  an  at¬ 
tempt  is  also  made  to  model  the  details  of  the  geologic 
environment  in  which  the  structure  is  to  exist.  For 
example,  the  in-si tu  stress  state,  the  strength  and  de- 
formability  of  the  different  rock  formations,  and  the 
frequency,  orientation,  and  strength  and  deformation 
characteristics  of  discontinuities  such  as  joints,  bedding 
planes,  and  faults  are  modeled  as  accurately  as  is  pos¬ 
sible  and  p r  u  tical.  The  structure  and  surrounding  geo¬ 
logic  environment  arc  envisioned  as  a  single  interacting 
unit  in  which  (lie  behavior  of  the  structure  itself  can¬ 
not  he  predicted  without  giving  due  consideration  to 
the  behavior  of  the  surrounding  geologic  environment. 

The  use  of  geomcchnnical  models  appears  to  he 
the  only  technique  available  analytical  or  experi¬ 
mental  for  determining  the  behavior  of  underground 
openings  through  all  stages  of  loading  and  deformation, 
both  elastic  and  inelastic,  up  to  failure. 

The  theoretical  basis  of  model  studies  in  general 


(loi  example,  Murp.iy  '1  Langluiar,s )  and  of  structural 
and  geomechanieal  models  in  particular  (for  example, 
Preeec  and  Davies.''  Rocha.7  Fumagalli."  Mandcl,**) 
has  been  well  established,  flic  basis  for  geomechanieal 
model  studies  of  underground  openings  in  rock  (not 
considering  tuae-deoendent  behavior)  lias  been  dis¬ 
cussed  and  developed  to  varying  degrees  by  experi¬ 
menters  sindi  as  Barron  and  Larocque,10  Fvcrling,1  1 
Hobbs.1*  lloek,in  and  Fumagalli.1 4 

Some  critical  aspects  of  similitude  requirements 
such  as  boundary  loading  conditions  and  model  materi¬ 
al  properties  generally  have  not  been  adequately  satis¬ 
fied  by  tlie  invcstigalois  above,  however,  and  their 

*  G.  Murphy,  Similitude  in  engineering g  (Ronald  Press,  1950). 

4  ILL.  Langhaar.  D.mensional  Analysis  and  Theory  of  Models 
(John  Wiley  «*t  Sons.  1951). 

‘  B.W.  Prcece  and  J.W.  Davies,  Models  for  Structural  Concrete 
(C.R.  Hooks  Ltd..  London.  1964). 

I  M.  Rocha,  "Model  Tests  in  Portugal.”  Civ.  lingr.  and  Pub. 
Work  Raw,  Vo!  53.  No.  619  (January  1958)  pp  49-53.  and 
No.  620  (February  (958)  (>p  179-182:  ‘'Structural  Model 
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work  is  of  limited  value  in  a  general  understanding  of 
tin.'  influence  of  tire  variables  affecting  the  behavior  of 
underground  openings  in  rock,  l  or  example,  in  the 
studies  cited  above  the  rock  tunnels  studied  were  sub¬ 
jected  to  plane  stress  rather  than  the  plane  strain  load¬ 
ing  which  tunnels  are  subjected  to  in  the  field.  In 
addition,  the  model  materials  used  hv  the  investigators 
listed  above  had  angles  ol  internal  friction  ranging  from 
5  to  20  under  significant  stiesses  and  thus  would  lead 
to  an  underestimation  ol  the  strength  of  rock  materials 
around  tunnels  in  the  field. 

Holier  and  Hendron* s  have  overcome  these  diffi¬ 
culties  hv  developing  a  device  to  subject  model  tunnels 
to  plane  strain  loading  and  by  developing  a  model 
material  with  an  angle  of  internal  friction  comparable 
to  many  real  rock  materials.  In  addition,  techniques 
have  been  developed  for  acquiring  detailed  quantitative 
measurements  of  the  behavior  of  the  rock  mass  behind 
the  tunnel  wall.16  Recent  developments  also  enable 
the  modeling  of  the  interaction  of  structural  linings 
with  jointed  rock  masses.' 1 

Scope  of  Study.  In  this  study  a  number  of  models 
were  tested  to  study  their  behavior.  The  models  were 
lined  with  plexiglass  liners  of  varying  thicknesses  to 
simulate  reinforced  concrete  tunnel  liners  in  a  jointed 
rock  mass.  The  effect  of  the  ratio  of  the  tunnel  diame¬ 
ter  to  joint  spacing  on  structural  behavior  was  also 
studied  by  testing  three  different  size  openings  in  the 
mode/  rock  mass  composed  on  I -in.  square  joint 
blocks.  The  structural  behavior  of  the  model  was  as¬ 
sessed  by  means  of  quantitative  measurements  such  as 
clip  gage  measurements  of  diameter  changes  of  the 
liner,  and  model  extensometers  to  measure  displace¬ 
ments  in  the  rock  mass  behind  the  tunnel  wall. 
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I  he  similitude  requirements  governing  the  choice 
ol  mode1  rock  material  and  the  model  structural  liner 
are  given  in  Appendix  A.  Although  this  analysis  was 
developed  in  previous  studies  by  Hcuer  and  Hendron,  it 
is  repeated  in  Appendix  A  to  give  the  reader  a  better 
understanding  of  the  basis  for  the  choice  of  the  struc¬ 
tural  linings  tested  in  this  study. 

The  methods  and  dev.uk  of  constructing  the 
jointed  models  are  described  in  Chapter  2,  which  also 
reviews  strength  properties  of  the  intact  model  material 
and  gives  strength  properties  of  the  joint  surfaces.  De¬ 
tails  of  the  "odel  instrumentation  arc  also  included. 
(Appendix  If  describes  the  model  loading  apparatus.) 
Chapter  3  discusses  selection  of  the  model  structural 
liners  and  icp  vis  results  of  the  experimental  measure¬ 
ments.  An  analysis  of  the  data  and  the  conclusions 
which  can  he  drawn  from  these  tests  is  in  Chapter  4. 


2  GEOMECHANICALMODELING 
TECHNIQUES 

The  model  material,  the  techniques  of  model  con¬ 
struction.  and  the  instrumentation  and  loading  equip¬ 
ment  used  on  this  contract  (DACA  23-70-C-0050)  were 
all  developed  on  a  previous  contract  (DACA 
39-(>7-C-0009)  with  the  U.S.  Army  Engineer  Water¬ 
ways  Experiment  Station.  Complete  discussions  of 
these  earlier  phases  of  the  project  arc  given  by  Hcuer 
and  Hendron.  Ileuer.  and  Hendron  cl  al.  The  remainder 
of  this  chapter  is  a  brief  discussion  of  the  modeling 
techniques  used  in  this  phase  of  the  project. 

Development  of  Jointed  Models.  The  loading  frame 
used  was  designed  to  test  24”  x  24"  x  8”  models  in 
plane  strain  (no  strain  along  the  axis  of  the  tunnel). 
The  model  tunnels  tested  on  this  study  were  4.  6,  and 
8  inches  in  diameter  and  were  drilled  through  the 
center  of  the  24”  x  24"  faces.  Thus,  the  model  tunnel 
simulates  a  section  in  the  axial  direction.  All  of  Ihc 
jointed  models  were  constructed  to  have  two  sets  of 
mutually  perpendicular  joints  oriented  parallel  to  the 
tunnel  axis.  Figures  1  and  2  show  the  joint  configura¬ 
tions  used  in  the  models  tested  in  this  study.  The 
models  were  tested  with  the  24"  x  24"  faces  hori¬ 
zontal.  and  thus  the  longitudinal  direction  is  vertical  in 
the  model  whereas  it  would  generally  be  the  horizontal 
direction  in  the  field.  The  models  were  tested  in  this 
orientation  because  it  greatly  simplified  the  design  of 
the  loading  apparatus. 


Figures  I  and  2  show  that  a  large  mmibcr  of  joint 
Mocks  were  required  for  Ihe  eonstruction  of  a  single 
model.  I  wo  possible  methods  could  be  used  to  manti- 
factuie  such  a  large  number  of  joint  blocks:  they  could 
either  be  cast  in  a  mold  to  the  proper  shape,  or  they 
could  be  sawed  'ut  of  larger  blocks  of  model  material. 
Bee  mse  ot  the  large  amount  of  time-consuming  work 
anticipated  in  a  sawing  process,  it  w;.s  dec.dcd  first  to 
try  molding  the  blocks  by  vibiating  a  sand-water* 
plaster  mix  in  a  mold.  The  anticipated  model  blocks 
would  be  required  to  have  a  l<  w  cohesion,  c,  so  they 
could  be  failed  by  the  testing  machine  and  a  high  angle 
>n  shearing  resistance,  0,  to  accurately  simulate  the 
properties  of  rock.  It  was  necessary  for  the  blocks  to 
have  a  very  dense  packing  of  sanu  grains  to  prevent 
collapse  of  their  structure  at  high  confining  pres¬ 
sures.  1  *  Also  it  was  desired  to  use  the  same  kind  of 
sand  and  plaster  in  the  vibrated  model  material  as  had 
been  used  in  the  compacted  model  material  used  in  the 
solid  model  blocks.1 0 

Attempts  to  make  joint  blocks  by  vibrating  materi¬ 
al  in  a  mold  proved  to  be  futile  because  the  blocks 
were  too  fragile  to  be  removed  from  the  mold.  (These 
were  2"  x  2"  x  8”  blocks).  The  failure  to  successfully 
extrude  the  vibrated  joint  blocks  was  due  largely  to  the 
very  low  cohesive  strength  of  the  material. 

After  attempts  at  molding  joint  blocks  failed  it 
was  decided  to  make  joint  blocks  by  sawing  them  out 
of  larger  compacted  blocks.  Steel  molds  20”  x  20”  x 
6”  were  used  to  compact  20”  x  20”  x  3”  blocks  using 
the  same  compaction  procedure  and  the  same  mix 
proportions  as  used  by  Hcucr  and  Hendron  on  24”  x 
24"  x  8”  solid  model  blocks  (Fig.  3).  A  decided  ad¬ 
vantage  of  this  procedure  is  that  the  intact  material  of 
ihe  joint  blocks  would  be  essentially  identical  to  the 
iniael  material  composing  the  solid  models  tested  pre¬ 
viously  by  Hcucr.  This  model  material  developed  by 
Hcucr  and  Hendron  is  probably  the  best  reported  to 
date  for  modeling  the  properties  of  rock. 

After  compaction,  the  blocks  arc  allowed  to  air 
dry  for  three  days,  then  put  into  an  oven  to  dry  at 
105°F  tor  about  a  week.  When  the  20”  x  20”  x  3” 
blocks  are  properly  cured,  the'-  arc  strong  enough  to  be 
handled  easily  without  breaking.  They  also  saw  very 
easily.  A  metal  surface  grinder  with  a  moving  table  was 

1 *  It. I  ..  lleiKT.  (U'onwchatiical  Model  Study. 
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converted  into  a  saw  for  accurately  cutting  joint  blocks 
(Fig.  4).  Diamond  blades  are  used  quite  successfully  for 
sawing  joint  blocks  with  this  machine  it  has  been  pos¬ 
sible  to  saw  blocks  as  small  as  1/2”  x  1/2”  x  8”.  A  jig 
was  also  made  to  lit  the  saw  for  cutting  the  triangular 
cross-section  blocks  used  around  the  edges  of  models 
which  have  joints  oriented  at  45°  to  the  principal  di¬ 
rections  of  loading  (Fig.  1 ). 

Since  exactly  the  same  material  was  used  in  the 
joint  blocks  as  was  used  by  Hcucr  and  Hendron  in  the 
solid  blocks,  a  new  series  of  material  properties  tests 
was  not  necessary.  The  standard  mix  is  made  in  the 
ratio  of  1.2/ I/O, '.01  ( water/plastcr/sand/rctardcr)  by 
weight.  The  plaster  used  is  white  molding  plaster.  The 
sand  is  the  line  fra*  lion  of  a  Pleistocene  sand  deposit 
obtained  from  the  Sangamon  River  valley  near  Ma¬ 
homet,  Illinois.  The  grain  size  distribution  of  the  Fine 
Sangamon  River  sand  is  shown  in  Fig.  5.  The  retarder 
used  is  sodium  phosphate  (Na2HP04)  in  the  dibasic 
anhydrous  powder  form. 

The  sand,  placer  and  retarder  arc  mixed  together 
dry  for  about  5  minutes  in  a  concrete  mixer.  The  water 
is  then  added  while  the  mixer  is  running  and  the  batch 
is  mixed  wet  for  about  5  minutes.  When  the  wet  mix  is 
homogeneous,  it  is  placed  in  the  mold  in  aboui  1/2 -in. 
thick  layers  and  compacted  with  a  pneumatic  tamper 
by  the  same  method  used  by  Hcucr  and  Hendron. 

The  intact  shear  strength  properties  of  the  model 
material  arc  shown  in  Fig.  6.  The  angle  of  internal 
friction  is  <p  =  33°  and  the  unconfined  compressive 
strength  is  q„  s  555  psi. 

The  Mohr  failure  envelope  for  the  intact  material 
in  Fig.  6  is  essentially  a  straight  line  up  to  confining 
pressures  as  high  as  1 000  psi.  This  is  in  marked  contrast 
to  the  behavior  of  most  previous  model  materials 
which  approach  0-0°  behavior  al  high  pressure.  Since 
a  high  frictional  shearing  resistance  is  one  of  the  most 
important  properties  of  jointed  rock  masses,  it  is  es¬ 
sential,  that  a  model  rock  material  have  high  frictional 
resistance. 

A  scries  of  three  direct  shear  tests  were  run  on  2” 
x  6”  sawed  joint  surfaces  of  the  model  material.  These 
tests  were  conducted  in  the  direct  shear  machine  in  the 
University  of  Illinois  rock  mechanics  laboratory.  Tests 
were  run  at  normal  stresses  of  50  psi,  150  psi.  and  400 
ps  .  The  measured  maximum  shear  strength  in  each 
case  respectively  was  33.3  psi,  97.5  psi,  and  230  psi. 
These  three  points  arc  plotted  in  Fig.  7,  which  shows 
that  the  effective  angle  of  shearing  resistance  on  the 
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Figure  1.  Jointed  model  configuration  used  in  Joint  Block  r*6,  #7  and  -#8. 


Figure  2.  Jointed  model  coni  imitation  used  in  Joint  Block  #0  anu  ^10. 
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joint  mi  ibices  decreases  from  33“  to  2(>“  with  increas¬ 
ing  normal  pressures.  All  three  direct  shear  specimens 
had  Hat-top  shear  strength  vs.  deformation  curve.:  for  a 
given  normal  pressure.  In  all  three  cases,  the  residual 
shear  strength  alter  3  cm  of  slip  along  the  joint  was 
essentially  the  same  as  the  peak  shear  strength.  These 
tests  indicate  that  a  value  of  the  angle  of  shearing  re¬ 
sistance  for  use  in  an  analysis  of  the  jointed  models 
should  he  slightly  lower  than  the  value  obtained  in  the 
triaxial  tests  of  intact  samples  shown  in  Fig.  6.  For  the 
theoretical  elasto-plastic  analysis,  the  appropriate  angle 
of  frictional  resistance  should  be  taken  as  the  angle  of 
frictional  resistance  along  the  joints,  not  the  angle  of 
internal  friction  derived  from  triaxial  tests  on  intact 
samples  of  the  model  material  (Hendron  and  Aiver.20) 

The  sawing  tolerance  on  the  hlov.  s  is  about  4  .01 
in.  This  means  that  in  a  model  with  ~  in.  joint  spacing, 
if  all  of  the  blocks  on  one  row  are  0.01  in.  too  thin  and 
on  the  next  row  they  are  all  0.01  in.  too  thick,  the 
maximum  offset  of  the  joints  across  the  model  would 
accumulate  to  as  much  as  0.24  in.,  which  is  intolerable. 
The  test  blocks  must  therefore  be  constructed  by  se¬ 
lecting  the  blocks  so  that  they  fit  together  to  make 
straight  joint  lines  with  minimum  offsets  in  both  di¬ 
rections  (Fig.  8).  The  blocks  are  constructed  on  edge 
on  a  table  and  then  moved  block  by  block  into  the 
testing  machine.  Each  of  the  external  faces  of  tsc  con¬ 
structed  block  is  flattened  by  grinding  and  is  tin. rough¬ 
ly  cleaned  of  dust  with  compressed  air  before  placing  it 
in  the  testing  machine. 

The  jointed  models  are  placed  in  the  testing  ma¬ 
chine  on  a  friction  reducing  sandwich  composed  of  a 
sheet  of  teflon  plastic  sandwiched  between  two  sheets 
of  4  mil  polyethylene  plastic  placed  directly  on  the 
base  plate  of  the  testing  machine.  The  plastic  sheets  arc 
used  to  reduce  friction  between  the  model  and  the  base 
plate  of  the  testing  machine. 

When  the  model  is  constructed  in  the  testing 
machine,  the  loading  elements  are  put  in  place  and  a 
small  seating  load  of  about  25  psi  is  applied  in  both  the 
horizontal  and  vertical  directions.  With  the  scaling  load 
held  constant,  the  4-in.  diameter  tunnel  is  cored  and 
cleaned  out  thoroughly  with  a  vacuum  cleaner.  The 
joints  intersecting  the  tunnel  are  then  sealed  with  a 
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Figure  8.  Jointed  model  ready  for  test. 


small  bead  of  silastic  caulking  compound  which  is  al¬ 
lowed  to  cure  for  two  days.  The  tunnel  wall  is  then 
painted  with  SR-4  strain  gage  cement  for  waterproof¬ 
ing  purposes  and  the  cement  is  allowed  to  cure  for  one 
day.  The  instrumented  tunnel  liner  is  then  installed  in 
the  tunnel  and  the  base  of  the  liner  is  scaled  with  si¬ 
lastic.  When  the  silastic  has  cured,  the  liner  is  grouted- 
in  with  a  liquid  grout  consisting  of  one  part  water  to 
one  part  sulfaset  rock  bolt  cement  by  weight.  The 
grout  is  cured  for  one  day  and  then  the  loading  head  is 
placed  on  top  of  the  model  using  two  layers  of  4  mil 
polyethylene  sheet  and  a  layer  of  plaster  to  get  close 
contact  between  the  model  and  the  testing  head.  This 
procedure  for  placement  of  the  loading  head  is  the 
same  as  that  used  by  Meucr. 

Free-Field  Strain  Measurement.  In  the  jointed  blocks, 
strain  gage  measurement  on  the  intact  blocks  are  mean¬ 
ingless  as  measure  of  the  free- field  strain  of  the  model 
due  to  closure  along  the  joints.  Thus  buried  extenso- 
meters  were  used  in  the  jointed  models  to  measure  the 
average  relative  displacement  of  two  points  across  the 
block.  Average  strains  of  the  block  were  obtained  by 
dividing  the  relative  displacement  between  the  points 
by  the  distance  between  the  points,  (<  =  AL/L).  Fig.  9 
shows  the  locations  of  extensometers  and  their  identifi¬ 
cation  code. 

The  buried  extensometers  are  simply  metal  rods 
grouted  with  epoxv  into  holes  bored  into  the  model  to 
the  specified  depth.  The  extensometer  holes  are  bored 
with  a  masonry  bit.  Elastic  tubing  is  used  to  contain 
the  epoxy  until  it  is  extruded  by  pushing  theextenso- 
meler  into  position  in  the  model.  Figure  10  shows  a 
series  of  extensometers  ready  to  be  filled  with  epoxy 
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Locations  of  cxtcnsometers  in  jointed  test  l)locks. 


Figure  10.  Extensometers  ready  for  installation  in  test 
block. 


and  installed.  The  extensometers  are  composed  of  a 
piece  of  1/4-in.  diameter  metal  rod  with  a  seating  hole 
drilled  in  the  end.  Two  pieces  of  plastic  tubing  are  used 
to  make  a  cup  to  contain  the  epoxy  on  one  end  of  the 
extensometer  and  another  piece  is  used  as  a  spacer  to 
center  the  lod  in  the  hole  near  the  face  of  the  model 
(Fig.  I  la).  A  typical  overall  installation  of  two  buried 
extensometers  is  shown  in  Fig.  lib.  The  epo.xv  is  very 
viscous  and  will  not  (low  out  when  the  extensometers 
arc  held  in  a  horizontal  position.  Figure  12  shows  an 
extensometer  as  it  was  exposed  after  a  test.  The  model 
must  be  under  a  seating  load  when  the  extensometer 
holes  arc  drilled  and  when  the  extensometers  are 
grouted  in  place. 

The  movements  of  the  extensometers  arc  meas¬ 
ured  with  beryllium-copper  clip  gages  like  those  used 
to  measure  diameter  changes  in  the  tunnels  of  the  pre¬ 
vious  tests.2'  Tnese  gages  consist  of  5 -in.  curved  strips 
of  beryllium-copper  with  seating  points  on  the  ends. 
Each  strip  is  gaged  with  four  strain  gages  wired  on  a 
four-arm  bridge.  These  clip  gages  can  be  accurately  cali¬ 
brated  with  a  standard  strain  indicator  and  they  are 
linear  over  a  range  of  about  1  in.  deflection.  Figure  I  lb 
is  a  detail  of  the  setup  used  to  measure  the  movements 
of  the  buried  extensometers. 

Each  of  the  clip  gages  is  calibrated  before  and  after 
each  test  while  wired  to  the  same  terminals  used  during 


the  test.  There  are  small  changes  in  the  calibration  of 
the  gages  from  test  to  test  and  continual  recalibration 
is  necessary  to  detect  changes  in  the  system  before 
miming  a  test. 


3  TEST  RESULTS 


Table  1  summarizes  the  model  tests  completed  in 
this  study.  The  five  tests  completed  are  designated  as 
JB  #6  through  JB  #10.  JB  #6  was  tested  at  an  N  value 
of  2/3  as  a  continuation  of  the  first  series  of  jointed 
models  (JB  #1  through  JB  #5)  tested  under  a  previous 
contract  (DACA  39-67-C-0009).22  JB  #7  through  JB 
#10  constitute  the  second  related  series  of  jointed 
models  tested.  All  the  models  in  this  second  series  were 
tested  at  a  principal  stress  ratio  N  =  oj,/ov  =  1  so  that 
the  results  could  be  compared  to  an  elasto-plastic  theo¬ 
ry  with  dilatancy  developed  by  Hendron  and  Aiyer. 
Comparisons  between  this  theory  and  the  model  test 
results  are  presented  in  Chapter  6. 


The  usual  type  of  liner  of  interest  for  practical 

protective  structures  problems  is  a  reinforced  concrete 

integral  liner  in  rock.  For  example,  consider  an  18-in. 

thick  reinforced  concrete  liner  in  a  20- ft  diameter 

opening.  The  circumferential  stiffness  of  the  liner  would 

i  .  i  v.m  3x  lCf  psi  X  18  in. 

be  approximately  Et/R  >  |01J-j2,n/fl  -  450.000 

psi  and  the  jointed  rock  mass  might  have  an  effective 
plane  strain  deformation  modulus  of  about  Er  =  1  —2 
X  1 0h  psi.  These  figures  would  give  a  linear  stiffness  to 

,  ..  r  Et/R  450,000  psi 

rock  mass  stitfness  ratio  of  —  =  - — - - -hr-  .  = 

Er  1-2  X  10  psi 

0.45  -  0.225.  In  these  model  tests,  three  sizes  of 
plexiglass  liners  were  used  as  shown  in  Table  1 .  The 
circumferential  stiffness  of  these  liners  varied  from 
16,700  psi  to  25,000  psi,  whereas  the  range  of  model 
stiffness.  Em,  for  models  JB  #7  through  JB  #10  was 
61,000  psi  to  64,300  psi  (Table  1).  The  ratio  of  the 
circumferential  stiffness  of  the  liner,  Et/R,  to  the  plane 
strain  stiffness  of  the  model,  Em,  varied  from  0.261  to 
0.410  (Table  1).  Thus  the  models  had  a  circumferential 
stiffness  to  rock  mass  stiffness  ratio  similar  to  that  of 
the  prototype  structure  mentioned  above.  The  ratio  of 
the  bending  stiffness,  El/R3,  to  the  circumferential 
stiffness,  Et/R,  determines  whether  a  liner  will  fail  by 
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I'i^nrr  1  la.  Detail  ni  l>iuie<l  I'x'iNM'im.’k'f. 


1‘ii’iirc  I  lb.  Dei  ail  <>l  extensometer  measurement  setup. 


Figure  12.  Extcnsometer  exposed  after  test. 


buckling  or  by  simple  circumferential  compression  fail¬ 
ure.  This  ratio  is  dependent  only  on  the  value  of  t/R. 
The  value  of  t/R  for  the  prototype  cited  above  is  0.15 
while  values  for  the  liners  used  in  the  models  arc 
0.0625  and  0.046.  Thus  the  model  liners  arc  more  like¬ 
ly  to  buckle  than  is  the  18-in.  thick  concrete  liner  in 
the  20-ft  diameter  example  prototype  structure.  In  the 
model  icsts,  none  of  the  plexiglass  liners  buckled  and 
two  of  them  failed  in  circumferential  compression  as  a 
prototype  concrete  liner  would.  Thus  the  selection  of 
plexiglass  liners  prevented  a  buckling  mode  of  failure 
which  had  previously  been  experienced  with  aluminum 
liners.23 

The  remainder  of  this  chapter  consists  of  a  system¬ 
atic  presentation  and  discussion  of  the  data  obtained  in 
the  testing  of  the  models  JB  #6  through  JB  #10. 

Presentation  of  Test  Results 

Joint  Block  # 6 .  # 6  was  identical  to  JB  #524  ex¬ 
cept  that  a  plexiglass  liner  was  used  in  JB  and  an 
aluminum  liner  was  used  in  JB  #5.  The  change  to  plexi¬ 
glass  was  made  because  the  bending  stiffness  fEI/R3  = 
4.5  psi)  of  the  aluminum  liners  was  too  low  and  the 
circumferential  stiffness  (Et/R  =  175,000  psi)  was  too 
high  to  model  reinforced  concrete  liners  in  rock.  As  a 
result  the  aluminum  liners  buckled  before  the  ultimate 
thrust  capacity  of  the  aluminum  section  was  de¬ 
veloped.  Also,  the  circumferential  stiffness  of  the 


aluminum  liner  was  about  4.5  times  as  great  as  the 
stiffness  of  the  jointed  model  in  the  vertical  direction, 
which  caused  the  liner  to  attract  excessive  load  by 
arching.  Thus  it  was  desirable  to  increase  the  bending 
stiffness  of  the  liner  in  JB  #6  to  reduce  the  chance  of 
buckling  failure.  It  was  also  desirable  to  decrease  the 
thrust  stiffness  of  the  liner  to  a  value  less  than  the 
modulus  of  the  surrounding  model  rock  mass  so  that 
the  test  woidd  more  accurately  model  real  prototype 
structures  in  jointed  rock  masses.  A  4-in.  diameter 
plexiglass  liner  I /8-in.  thick  has  a  calculated  thrust 
stiffness  of  25,000  psi  and  a  calculated  bending  stiff¬ 
ness  of  8.15  psi.  These  stiffnesses  arc  based  upon  a 
Young’s  modulus  of  the  plexiglass,  E,  of  400,000  psi. 
The  value  of  E  was  determined  by  testing  an  8-in.  long 
sample  of  the  4-in.  diameter  liner  in  uniaxial  compres¬ 
sion  using  SR-4  strain  gages  to  measure  the  strains.  The 
average  vertical  model  stiffness  of  JB  #5  was  E,„  = 
38.000  psi  and  thus  the  thrust  stiffness  of  this  liner  was 
lower  than  the  model  stiffness  as  desired. 

With  the  aid  of  the  information  above  JB  #6  was 
constructed  using  1 -in.  joint  blocks  oriented  at  45°  to 
the  loading  directions;  and  a  4-in.  diameter  plexiglass 
liner  1  /8-in.  thick.  JB  #6  was  tested  at  a  principal  stress 
ratio  N  =  0|,/ov  =  2/3  to  compare  with  six  previous 
tests25  conducted  at  N  =  2/3. 

Figures  13  through  15  arc  summary  plots  of  the 
data  obtained  in  the  testing  of  JB  #6.  Figure  13  is  a 
plot  of  the  vertical  stress-strain  curves  of  the  model. 
The  stress-strain  curve  designated  as  “vertical-shallow” 
was  obtained  from  a  pair  of  extensometers  buried  3  in. 
deep  in  opposite  vertical  faces  of  the  model  as  shown 
in  Fig.  9  and  represents  the  average  strain  over  the 
central  18  in.  of  the  vertical  model  centerline.  Thus  the 
“vertical-deep”  stress-strain  curve  is  more  affected  by 
strain  concentrations  due  to  the  tunnel  and  the  “verti¬ 
cal-shallow”  stress-strain  curve  is  more  representative 
of  the  frec-field  strain  of  the  model.  After  initial  seat¬ 
ing  movements,  both  of  the  vertical  stress-strain  curves 
arc  almost  linear  up  to  a  vertical  model  stress  of  1200 
psi.  Above  1 200  psi  they  both  show  the  strain  to  in¬ 
crease  at  a  slightly  increasing  rate  with  pressure.  The 
separation  of  the  two  vertical  stress-strain  curves  indi¬ 
cates  that  the  movements  of  the  deep  extensometers 
were  affected  by  the  presence  of  the  tunnel. 

The  uamctrical  extensometers  in  JB  #6  were 
placed  at  angles  of  10°,  35°,  55°,  and  80°  from  the 
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Table  I 

Summary  of  Model  Blocks  Tested 


Test 

Block 

Tunnel 

Diameter 

(in.) 

Liner 

Properties* 

Joint 

Spacing 

(in.) 

Joint 

Orientation 

«h 

N  =  Fm(psi) 

Et/r 

Jll  »6 

4 

plexiglass 
t  -  0.1  25  in. 

Ilt/R  -  25.000  psi 
El/R3  =  8.15  psi 

1 

45° 

2/3  54,000 

0.463 

JB  #7 

4 

plexiglass 

1  =  0.125  in. 

Ilt/R  =  25,000  psi 
l'l/'R3  =  8.15  psi 

1 

45° 

1  64,300 

0.389 

JB  ttX 

6 

plexiglass 
t  =  0.125  in. 
l.t/R  -  16.700  psi 
El/R3  =  2.4  psi 

1 

45" 

1  64,000 

0.261 

JB  n 9 

6 

plexiglass 
t  =  0.125  in. 
l  l/R  =  16.700  psi 
111/ R 3  =  2.4  pst 

1 

90'“ 

1  62,500 

0.268 

JB*I1» 

K 

plexiglass 
t  =  0.250  it.. 
tll/R  =  25.000  psi 
l  l/R3  =  8.15  psi 

1 

00“ 

1  6 1 ,000 

0.410 

*  liner  stiffnesses  calculated  for  a  I  in.  length  of  the  liner. 


crown  and  invert  diameter  instead  of  the  usual  0°,  45°, 
and  90°  used  in  the  other  tests.  The  extensometers 
were  oriented  this  way  to  help  determine  the  actual 
shape  of  the  tunnel  liner  after  deformation.  Fig.  14 
is  a  plot  of  these  diametrical  strains  labeled  as  to  the 
angle  they  make  with  the  crown  and  invert  direction. 
These  curves  show  that  the  liner  was  actually  more 
flexible  than  the  model  as  even  the  springlinc  diameter 
was  compressed  in  this  test. 

Fig.  15  is  a  dimensionless  plot  of  the  diametrical 
strain  of  the  liner  on  four  different  diameters  as  a  func¬ 
tion  of  the  vertical  frcc-ficld  strain.  In  this  plot,  the 
vertical-shallow  stress-strain  curve  was  used  as  the  frec- 
ficld  stress-strain  curve.  Presentation  of  the  data  in  this 
manner  is  helpful  for  an  immediate  comparison  of  the 
diametrical  strain  AD/D  with  the  frcc-ficld  strain  in  the 
rock  mass.  Note  that  at  strains  less  then  1%  the  dia¬ 
metrical  strains  along  the  10°  diameter  and  the  35° 
diameter  arc  nearly  equal  to  the  frcc-ficld  strain.  As  the 
free  field  strain  increases  above  1%,  inelastic  action  of 
the  model  rock  mass  surrounding  the  liner  results  in 
the  diametrical  strains  increasing  at  an  increasing  rate 


with  respect  to  the  frcc-ficld  strain.  At  Tree-field  strains 
of  2  -3%  the  diametrical  strains  along  these  two  diame¬ 
ters  arc  about  1.5— 1.8  as  great  as  the  frcc-ficld  strain. 
The  diametrical  strains  on  the  55°  diameter  and  the 
80°  diameter  are  less  than  the  frec-ficld  strain  because 
of  ovaling.  At  higher  loads,  the  80°  diameter  would 
have  shown  an  actual  increase  over  its  original  length 
due  to  ovaling  of  the  liner  caused  by  testing  at  N  =  2/3. 

Joint  Block  # 7 .  The  second  series  of  jointed 
models  (JB  #7  through  JB  #10)  was  tested  at  a  princi¬ 
pal  stress  ratio  N  =  < q,/ov  =  1  specifically  to  compare 
the  test  results  to  a  theoretic  ll  elasto-plastic  analysis 
with  dilatancy  which  was  developed  by  Hendron  and 
Aiycr.  The  theory  is  valid  for  calculating  stresses  and 
strains  around  a  cylindrical  opening  in  a  homogeneous 
daslo- plastic  material  with  dilatancy  under  symmetri¬ 
cal  plane-strain  conditions  for  a  principal  stress  ratio  N 
=  oh/ov  =  1.0.  The  required  loading  conditions  arc  sat¬ 
isfied  in  (his  scries  of  model  tests.  One  purpose  of  these 
model  tests  is  to  provide  test  data  lor  tunnels  in  jointed 
media  so  that  appropriate  values  of  shear  strength  pa¬ 
rameters  for  jointed  ma.scs  may  be  determined  for  use 
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in  litis  theory  which  will  yield  realistic  results  lor  tun¬ 
nels  in  jointed  masses.  This  theory  has  already  been 
used  to  predict  the  non-linear  diameter  changes  of 
tunnel  liners  in  solid  models  of  the  same  model  materi¬ 
al  where  the  models  were  subjected  to  stresses  high 
enough  to  cause  inelastic  behavior  around  tlte  tun¬ 
nel.3'' 


JB  37  was  constructed  of  1  in.  joint  blocks  ori¬ 
ented  at  45°  to  the  loading  directions  The  tunnel  was 
4  in.  diameter  and  was  lined  with  a  l/s-in.  thick  plexi¬ 
glass  liner.  The  circumferential  stiffness,  Et/R.  of  this 
liner  was  25,000  psi  and  the  bending  stiffness,  El/F3, 
was  <S.  1 5  psi.  JB  #7  was  tested  at  a  principal  stress  ratio 
N  =  0|,/ov  =  1  to  a  peak  model  pressure  of  1420  psi. 
Figures  16  and  17  arc  summary  plots  of  the  data  ob¬ 
tained  in  the  testing  of  JB  #7. 

Fig.  16  is  a  plot  of  the  average  model  stress- 
strain  curves  of  JB  #7  as  measured  by  four  separate 
pairs  of  external  extensometers.  The  shallow  extenso- 
meters  measured  over  a  gage  length  of  18  in.  while  the 
deep  extensometers  measured  over  a  gage  length  of  14 
in.  An  average  of  these  four  stress-strain  curves  of  JB 
#7  gives  an  average  plane-strain  deformation  modulus, 
Em,  of  the  jointed  mass  of  64.300  psi.  Thus  the  ratio 
of  the  thrust  stiffness  of  the  liner  to  the  stiffness  of  the 
Et/R 

model  — r —  was  0.388.  The  vertical  deformation 
*-ni 

modulus  of  JB  #6  was  54.000  psi.  JB  #7  was  stiffer 
than  JB  #6  because  JB  #7  was  tested  at  N  -  1  while  JB 
36  was  tested  at  N  =  2/3 .  Under  testing  at  N  =  1,  the 
vertical  plane-strain  stiffness  of  identical  models  is 
higher  due  both  to  increased  confinement  of  the  model 
and  to  a  lack  of  shearing  deformations  along  the  joint 
surfaces.  If  these  curves  (Fig.  16)  were  all  true  free¬ 
hold  stress-strain  curves  (unaffected  bv  the  proscncc  of 
the  liner),  then  they  should  all  be  coincident  for  testing 
at  N  =  1.  In  Fig.  16  the  deep  stress-strain  curves  show 
more  strain  than  the  corresponding  shallow  stress-strain 
curves  because  of  the  effects  of  the  tunnel.  Also  in  Fig. 
16  the  horizontal  stress-strain  curves  show  about  15% 
more  strain  than  the  vertical  stress-strain  curves.  This 
relationship  was  not  expected  and  can  be  explained 
only  by  experimental  error. 

Fig.  17  is  a  plot  of  the  diametrical  strain  of  the 
liner  as  a  function  of  model  pressure  at  four  different 
diameters  spaced  45°  apart  in  the  liner.  These  four 
curves  arc  very  nearly  coincident  thus  indicaiing  that 
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the  model  was  behaving  almost  like  a  homogeneous 
isotropic  medium  in  spite  of  the  joints.  Note  that  up  to 
a  pressure  of  about  600  psi  the  behavior  of  the  liner 
was  essentially  clastic  and  at  pressures  above  600  psi 
the  liner  strains  began  to  increase  at  an  increasing  rate 
with  pressure.  This  i  idicatcs  inelastic  movement  and 
yielding  of  the  rock  mass  around  the  tunnel  liner. 

Joint  Block  #6’.  JB  #3  was  constructed  with  a  1- 
in.  joint  spacing  with  the  joints  oriented  at  45°  to  the 
loading  directions.  The  tunnel  in  JB  #8  was  6  in.  in 
diameter  and  was  lined  with  a  1/8-in.  thick  plexiglass 
liner.  The  calculated  circumferential  stiffness,  Et/R,  of 
this  liner  was  16,700  psi  and  the  calculated  bending 
stiffness,  EI/R3,  was  2.4  psi.  JB  #8  was  tested  at  a 
principal  stiess  ratio  N  =  0|,/ov  =  1  to  a  peak  model 
pressure  of  about  1320  psi.  At  this  pressure  the  liner 
began  to  collapse  and  the  test  was  terminated  to  pre¬ 
vent  destruction  of  the  clip  gages  inside  the  tunnel 
liner.  Figures  18  and  19  arc  summary  plots  of  the  data 
obtained  in  the  testing  of  JB  #8. 


Fig.  IS  is  a  plot  of  the  average  model  stress- 
strain  curves  of  JB  #8  as  measured  by  four  separate 
oairs  of  external  extensometers.  The  shallow  extenso- 
nicters  measured  over  a  gage  length  of  18  in.  while  the 
deep  extensometers  measured  over  a  gage  length  of  14 
in.  An  average  of  these  four  stress-strain  curves  of  JB 
#8  gives  a  deformation  modulus  of  the  jointed  mass, 
Em,  of  64,000  psi.  Thus  the  ratio  of  the  calculated 
circumferential  stiffness  of  the  liner  to  the  stiffness  of 


the  model 


Et/R 


for  JB  #8  was  9.260  as  compared  to 


0.388  for  JB  #7.  Fig.  18  shows  that  again  in  JB  #8 
the  horizontal  strains  of  the  model  were  slightly  greater 
than  the  vertical  strains  This  difference  must  again  be 
explained  by  experimental  error;  possibly  a  systematic 
error  since  it  is  in  the  same  direction  as  it  was  in  JB  #7. 
In  JB  *8,  the  deep  extensometers  again  measured  more 
strain  than  the  shallow  extensometers.  This  strain  gra¬ 
dient  is  due  to  the  strain  concentration  caused  by  the 
tunnel. 


Fig.  19  is  a  plot  of  the  diametrical  strain  of  the 
liner  as  a  function  of  model  pressure  at  four  different 
diameters  spaced  45°  apart  in  the  liver.  These  four 
curves  arc  nearly  coincident  thus  indicating  that  the 
model  was  behaving  almost  like  a  homogeneous  iso¬ 
tropic  medium  in  spite  of  the  joints.  The  liner  sus¬ 
tained  the  model  pressure  of  1270  psi  continuously 
without  further  significant  deformation  in  spite  of  the 
fact  that  it  was  on  the  verge  of  failure  and  actually 
failed  al  1 320  psi. 
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Figure  16.  Stress-strain  curves  for  JB  #7. 
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Model  Pressure 
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Figure  17.  Diametrical  strain  of  tunnel  liner  at  lour  different  diameters  for  JB  #7. 
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Figure  18.  Stress-strain  curves  for  JB  #8. 
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Model  Pressure  o- =  cr  (psi) 
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Figure  19.  Diametrical  strain  of  tunnel  liner  at  four  different  diameters  for  JB  #8. 
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Joint  Block  #9.  JB  #9  was  constructed  with  a  I- 
in.  lomt  spacing  with  the  joints  mien  ted  parallel  to  the 
loading  directions.  The  tunnel  in  JB  #9  was  6  in.  in 
diameter  and  was  lined  with  a  PH  in.  thick  plexiglass 
huer.  I  he  calculated  circumferential  stiffness  of  this 
liner,  l-.l/R,  was  !f>,700  psi  and  the  calculated  bending 
stillness.  Fl/R\  was  2.4  psi.  JB  #<>  was  tested  at  a 
principal  stress  ratio  N  =  oj,/r jv  =  |  in  a  p.  ;ak  model 
pressure  ot  about  1410  psi.  at  which  piessure  the  liner 
suddenly  cracked  loudly  and  the  test  was  terminated  to 
prevent  complete  collapse  of  the  tunnel  and  conse¬ 
quent  destruction  of  the  clip  gages  inside  the  tunnel 
liner.  Figures  20  and  21  are  summary  plots  of  the  data 
obtained  in  the  testing  of  JB  #0. 


Fig.  20  is  a  plot  of  tire  average  model  stress- 
strain  curves  o!  JB  =9  as  measured  by  four  separate 
pairs  ol  external  extensometers.  The  shallow  extenso- 
meters  measured  over  a  gage  length  of  IX  in.  while  the 
deep  extensometers  measured  over  a  gage  length  of  14 
in.  An  average  of  these  four  stress-strain  curves  of  JB 
#0  gives  an  average  plane-strain  deformation  modulus 
of  the  jointed  mass  of  about  62,500  psi.  Thus  the  ratio 
of  the  calculated  circumferential  stiffness  of  the  liner 

to  the  stiffness  of  the  model,  -il— for  JB  =9  was 
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0.267  as  compared  to  0.260  for  JB  =8.  Fig.  20 
shows  that  the  horizontal  strains  in  JB  =9  were  again 
slightly  greater  than  the  vertical  strains.  This  difference 
is  almost  surely  a  systematic  experimental  error.  The 
loading  frame  and  the  hydraulic  elements  of  the  load¬ 
ing  system  have  been  checked  and  no  cause  for  such  a 
systematic  error  was  found.  In  JB  #9  the  deep  extenso- 
meters  again  measured  more  strain  than  the  shallow 
extensometers  due  to  the  strain  gradients  near  the 
tunnel. 


Fig.  21  is  a  plot  of  the  diametrical  strain  of  the 
tunnel  liner  as  a  function  of  model  pressure  on  three 
different  diameters.  The  change  in  slope  of  each  curve 
at  a  stress  level  of  400  psi  indicates  that  inelastic  action 
in  the  rock  mass  surrounding  the  tunnel  was  initiated 
at  this  stress  level.  Inelastic  action  also  started  at  about 
400  psi  in  JB  #8.  Inelastic  action  began  at  about 
600  700  psi  in  JB  #7  which  had  a  4-in.  tunnel.  Since 
models  JB  #8  and  JB  *9  had  6-in.  diameter  tunnels 
and  yielded  at  lower  stress  levels  than  JB  -7.  the  de¬ 
crease  in  rock  mass  strength  with  increasing  tunnel  size 
is  illustrated  by  these  tests.  This  point  will  be  discussed 
in  the  next  chapter. 

Joint  Block  # 10 .  JB  #10  was  constructed  with  a 
I  in.  joint  spacing  with  the  joints  oriented  parallel  to 


the  loading  directions,  flic  tunnel  in  JB  =10  was  h  m 
in  diameter  and  was  lined  with  a  l/4-ni.  thick  plexiglass 
liner.  Using  a  Young's  modulus  of  the  plexiglass  in  the 
circumferential  direction  ol  400,000  psi,  the  calculated 
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circumferential  stillness  of  the  liner, _i— — .'.was  25.000 

Fm 

psi  and  the  calculated  bending  stiffness.  Fl/R’.  was 
K  I  5  psi.  JB  #10  was  tested  at  a  principal  stress  ratio  N 
=  i>ii/'\  I  to  a  peat-  model  pressure  of  14  20  pm 
without  failing.  Figures  22  and  24  are  summary  plots 
of  the  data  obtained  in  the  testing  of  JB  =10. 

Fig.  22  is  a  plot  of  the  average  model  stress- 
strain  curves  of  JB  =10  as  measured  by  four  pairs  of 
external  extensometers.  An  average  of  these  four  siress- 
strain  curves  of  JB  =  10  gives  an  average  plane-strain 
deformation  modulus  of  the  jointed  mass  oi  about 
61.000  psi.  Thus  the  ratio  of  the  calculated  thrust  still¬ 
ness  of  the  tunnel  liner  to  the  aveiage  stiffness  of  the 
model,  Ft/R  was  0.410.  Fig.  22  shows  that  again  in 
JB  #10  the  horizontal  strains  were  greatci  titan  the 
vertical  strains.  The  deep  extensometers  also  measured 
more  strain  than  the  shallow  extensometers  because  of 
the  strain  concentrations  near  the  tunnel. 

Fig.  25  is  a  plot  of  tlte  diametrical  strain  of  the 
tunnel  liner  as  a  function  of  the  model  pressure  along 
four  different  tunnel  diameters.  The  data  shown  in  this 
figure  is  different  than  for  the  previous  joint  blocks 
discussed  because  there  was  no  definite  pressure  at 
which  yielding  was  apparent.  Furthermore  the  average 
diametrical  strain  at  1400  psi  was  about  3%  for  this  8- 
in.  tunnel  liner,  whereas  the  4-in.  tunnel  tested  in  JB 
#7  showed  about  5r/e  strain  al  comparable  slress  levels. 
This  was  not  expected  because  both  tunnel  liners  were 
selected  for  having  the  same  circumferential  stiffness. 
Et/R,  and  the  8-in.  tunnel  liner  should  have  shown 
larger  strains  because  the  ratio  ol'  tunnel  diameter  to 
joint  spacing  was  8.0  rather  than  the  value  of  4.0  used 
in  JB  =  7. 

Comparisons  of  Test  Data  From  JB  #7  through  JB 
#10.  Figures  24  through  26  arc  summary  plots  of 
average  data  obtained  in  the  tests  of  JB  #7  through  JB 
#10.  The  relative  behavior  of  the  various  tunnel  liners 
may  be  compared  directly  on  these  figures.  All  of  these 
models  (JB  =7  through  JB  =10)  were  tested  al  a  princi¬ 
pal  stress  ratio  N  =  <J|,/ov  =  I  and  all  had  joints  spaced 
at  1  in.  in  two  mutually  perpendicular  directions.  Thus 
all  variables  in  these  four  tests  were  held  constant  ex¬ 
cept  tlte  tunnel  size  and  linei  stiffness. 

Fig.  24  is  a  pint  of  the  average  stress-strain  curve 
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Figure  21.  Diametrical  strain  of  tunnel  liner  :>t  three  different  diameters  for  !.M 
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Figure  23.  Diametrical  strain  of  tunnel  liner  at  four  different  diameters  for  JB  ^10. 


lor  each  model  obtained  by  averaging  ail  four  of  the 
measured  stress-strain  curves.  Each  curve  is  adjusted  to 
zero  strain  at  a  model  stress  of  125  psi  to  eliminate  the 
elfects  of  nonuniform  seating  movements.  If  the  four 
models  were  all  properly  constructed  with  the  same 
tolerances  and  model  properties,  then  for  a  loading  at  a 
principal  stress  ratio  N  =  0|,/ov  =  I  they  should  all  have 
the  same  stress-strain  curve.  The  four  curves  on  Fig.  24 
are  remarkably  similar;  each  shows  additional  sealing 
movements  up  to  a  model  pressure  of  about  300  psi 
and  then  is  essentially  linear  with  a  slope  ranging  from 
61,000  psi  to  64,300  psi.  The  two  models  with  joints 
oriented  at  00°  to  the  loading  directions  (JB  *9  and  JB 
#10)  show  slightly  greater  strains  at  similar  model  pres¬ 
sures  than  the  two  with  joints  oriented  at  45  ’. 

Fig.  25  is  a  plot  of  the  average  diametrical  strain 
of  the  tunnel  liners  as  a  function  ol  model  pressure  for 
each  of  the  four  models  (JB  #7  through  JB  -10).  All 
of  these  models  were  tested  at  N  =  I.  and  all  had  very 
nearly  similar  average  stress-strain  curves  as  seen  in  Fig. 
24.  Thus,  the  major  differences  between  the  four 
curves  on  Fig.  25  should  be  caused  by  differences  in 
the  stiffness  of  the  liners  and  differences  in  the  ratio  of 
joint  spacing  to  tunnel  diameter.  The  joint  spacing  was 
1  in.  in  all  four  cases,  so  the  ratio  of  joint  spacing,  s.  to 
tunnel  diameter.  D.  was  1/4  for  JB  #7,  i  /6  for  JB  #8 
and  JB  #9,  and  1/8  for  JB  #10.  It  has  been  suggested 
by  Hcndron  and  Aiycr  on  the  basis  of  field  observa¬ 
tions  that  a  decreasing  value  of  s/D  should  cause  in¬ 
creased  load  on  a  tunnel  liner  at  the  same  value  of  free 
field  stress. 

In  Fig.  25,  the  curves  for  JB  #8  and  JB  £9  are 
almost  exactly  coincident  up  to  a  model  pressure  of 
1000  psi.  These  two  models  had  the  same  stiffness 
(Et/R  =  io,700  psi)  and  diameter  liner  and  differed 
only  in  joint  orientation.  Joint  orientation  should  have 
no  effect  for  testing  at  N  -  0|,/ov  =  1  if  the  boundary 
conditions  in  the  testing  machine  arc  as  assumed. 

JB  #7  had  a  4-in.  diameter  liner  t/8-in.  thick 
(Et/R  »  25,000  psi).  Thus  JB  *7  and  JB  *10  both  had 
the  same  calculated  stiffness  liner  and  about  the  same 
model  stiffness.  This  is  evidenced  by  the  fact  that  the 
plot  of  pressure  versus  diametrical  strain  lor  these  tests 
arc  nearly  identical  straight  lines  in  Fig.  25  up  to  a 
model  pressure  level  of  550  psi.  Note  also  that  the  test 
results  from  JB  #8  and  JB  #9  (Fig.  25)  arc  also  identi¬ 
cal  straight  lines  up  to  about  500  ps'  at  a  lower  slope 
than  tests  JB  #7  and  JB  #10  because  the  liner  stiffness 
in  JB  #8  and  JB  #9  was  onlv  16.700  psi.  Above  a 


pressure  level  of  550  psi  the  curves  for  JB  #7  and  JB 
#10  (Fig.  25)  diverge  as  the  behavior  of  the  jointed 
mass  surrounding  the  tunnel  becomes  inelastic  but  the 
trend  was  surprising  because  it  was  thought  that  the  8 
in.  tunnel  (JB  #10)  would  show  larger  diametrical 
strains  than  the  4  in.  Hum'.-'  (JB  #7)  with  the  same 
liner  stiffness. 

In  Fig.  26  the  curves  of  diametrical  strain  of  the 
tunnel  liner  versus  the  average  frcc-ficld  strain  for  joint 
blocks  7  through  10  am  shown.  Note  that  at  free  field 
strains  below  0.8%  the  strain  concentration  factor 

(AD/O)  tjic  |UllIlL.i  |in,.rs  in  JB  #7  and  JB  #10 
cff 

(Ei/R  =  25.000  psi )  are  lower  than  the  strain  concen¬ 
tration  factors  for  the  tunnel  liners  tested  in  JB  #8  and 
JB  #9  (Et/R  =  16,700  psi).  At  higher  values  of  the 
frcc-ficld  strain  the  lest  results  of  JB  #7,  #8.  and  #9 
are  nearly  identical  in  having  strain  concentration  fac¬ 
tors  of  about  2  and  2.5  at  free  field  strains  of 
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1.2%  and  2%  respectively.  Test  block  JB  *10  did  not 
plot  with  the  other  icsnlts  because  the  true  frcc-ficld 
strains  of  the  model  were  not  measured  in  JB  *10.  The 
8-in.  diameter  tunnel  was  so  large  with  icspecl  to  the 
size  of  the  model  that  die  measured  strains  were  af¬ 
fected  by  the  large  diametrical  strains  of  the  tunnel. 


4  ANALYSIS  OF  TEST  RESULTS 

Hcndron  and  Aiycr  have  analyzed  the  stresses  and 
.-.trains  around  a  cylindrical  tunnel  in  a n  elasto-plastic 
material  with  dilatancy  for  a  uniform  stress  field  (N  = 
oi, /«v  =  I  )  under  plane  strain  conditions.  The  loading 
conditions  assumed  in  the  analysis  cited  above  are  ex¬ 
actly  the  same  as  the  conditions  in  the  model  tests 
conducted  on  this  study.  The  main  problem  in  apply¬ 
ing  the  theory  to  the  jointed  models  is  to  determine 
what  strength  properties  should  be  used  in  the  analysis 
to  represent  a  jointed  mass.  Ilendron  and  Aiycr  give  a 
correlation  between  the  ratio  o„/q„  and  the  ratio  D/s 
lot  actual  tunnels  in  jointed  rock  masses  (Fig.  27). 
where  <r„  is  the  equivalent  unconfincd  compressive 
strength  ot  the  jointed  rock  mass  around  the  tunnel.  q„ 
is  the  intact  unconfined  compressive  strength  of  the 
rock,  !)  is  the  tunnel  diameter  and  s  is  the  effective 
joint  spacing.  The  correlation  is  from  field  measure¬ 
ments  of  tunnel  behavior.  The  models  tested  under  this 
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Figure  25.  Average  diametrical  strain  of  tunnel  liner  as  a  function  of  model  pressure  for  JB  -7,  JB  #8,  JB  #9  and 
JB  #10. 
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Figure  26.  Dimensionless  plot  of  average  diametrical  strain  of  tunnel  liner  as  a  function  of  average  vertical  strain  of 
model  for  JB  #7,  JB  #8,  JB  #9  and  JB  *10. 
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Figure  27.  Ratio  of  insitu  strength  to  laboratory  strength  as  a  function  of  tunnel  diameter  to  joint  spacing. 


contract  had  values  of  D/s  of  4,  6,  and  8  with  respec¬ 
tive  values  of  ou/qu  obtained  from  Fig.  27  of  0.65,  0.5. 
and  0.35.  Since  the  model  material  has  an  intact  un¬ 
confined  compressive  strength  of  550  psi,  these  values 
yield  ou  values  of  360  psi,  275  psi.  and  190  psi  respec¬ 
tively.  The  field  correlations  were  from  real  tunnels 
where  the  joints  were  not  exactly  regular  plane  surfaces 
nor  were  they  spaetd  evenly,  and  all  oriented  in  two 
sets  exactly  parallel  to  the  tunnel  axis.  Thus  the  field 
situation  is  not  in  general  as  unfavorable  as  the  situ¬ 
ation  in  the  models  and  a  lower  value  of  ou/qu  would 
be  expected  for  the  models  than  for  the  field  con¬ 
ditions. 

From  the  model  tests  conducted  on  this  study 
values  of  ou  can  be  calculated  using  the  theory  above 
in  conjunction  with:  (1)  measured  values  of  AD/D  = 
e0,  ( 2)  measured  values  of  Em,  (3)  calculated  values  of 
the  radial  stresses  between  the  liner  and  the  medium 
from  the  liner  properties  and  measured  values  of  AD/D 


=  e0  and  (4)  known  or  estimated  values  of  u  and  <t>  for 
the  jointed  mass.  Items  (11(2)  and  (3)  above  are  easily 
obtained  from  measurements  on  the  model  during  test¬ 
ing  but  t he  values  of  v  and  <p  for  the  jointed  mass  arc 
not  known  precisely.  The  value  of  Poisson's  ratio  can 
be  taken  as  about  0.25  based  on  previous  experience 
but  the  estimation  of  this  value  is  not  critical  since  the 
results  of  the  theory  cited  above  arc  not  very  sensitive 
to  the  selected  value  of  v.  The  value  of  the  angle  of 
shearing  resistance  for  the  rock  mass  may  be  conserva¬ 
tively  taken  as  the  angle  of  friction  along  the  joint 
surfaces,  7 

The  value  of  the  angle  of  shearing  resistance  along 
the  joints,  0j,  is  estimated  from  the  results  of  direct 
shear  tests  on  sawed  joint  surfaces  shown  in  Fig.  7.  The 
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results  ot  the  direct  shear  tests  on  joint  surfaces  show 
that  the  value  ol  0j  decreases  as  the  normal  pressure  on 
the  joint  surface  increases.  Thus  for  plastic  behavior  of 
the  jointed  model  we  should  expect  an  effective  value 
of  less  than  the  angle  of  internal  friction,  which  was 
measured  as  33°  by  triaxial  tests. 

The  parameters  used  in  the  Hendron-Aiycr  elasto- 
plastic  analysis  are  F.m,  /u,  </>j,  ou,  eH.  rr,  and  p„.  Where 
Em,  n,  and  ou  are  effective  plane  strain  properties  of 
the  jointed  model  mass,  e0  and  of  are  strain  and  stress 
conditions  at  the  point  of  interest  and  p„  is  the  free 
field  stress  in  the  model  (po  =  ctj,  ~  av  for  the  model 
tests  reported  herein).  Theoretical  curves  of  the  dia¬ 
metrical  strain  of  the  tunnel  liner,  AD/D.  as  a  function 
of  the  model  pressure.  oh  =  ov  =  p„.  have  been  calcu¬ 
lated  using  measured  values  of  Em  and  assumed  values 
of  /i,  and  ou.  With  known  values  of  the  liner  proper¬ 
ties,  Fig.  28.  or  can  be  calculated  as  a  function  of  e0  = 
AD/D  of  the  liner.  Then  a  value  of  p<(  can  be  calculated 
from  the  theory  and  theoretical  curves  of  AD/D  =  c0 
versus  p„  can  be  plotted.  The  actual  experimental 


curves  from  the  model  tests  can  be  plotted  on  the  same 
plot  to  get  a  comparison  between  theory  and  experi¬ 
ment.  Such  curves  are  plotted  in  Figures  29  and  30. 

JO  #S  and  JU  #0  were  identical  in  all  respects 
except  that  they  had  different  joint  orientations.  Joint 
orientation  should  have  no  effect  for  testing  at  a  princi¬ 
pal  stress  ratio  N  =  </| ,/ov  =  1.  Fig.  24  shows  that 
indeed  JB  ~8  and  JB  #9  did  behave  almost  identically 
on  a  plot  of  diametrical  strain  of  tunnel  liner  AD/D  as 
a  function  of  model  pressuri  oj,  =  av  =  p„.  In  Fig.  29, 
the  results  of  the  model  tests  JB  #8  and  JB  #9  are 
plotted  together  with  theoretical  curves  calculated 
from  the  Hendron-Aiyer  analysis  assuming  the  liner 
properties  shown  on  Fig.  28,  and  the  model  properties 
listed  on  Fig.  29.  The  curves  for  JB  #8  and  JB  #9  are 
essentially  identical  up  to  a  model  pressure  of  1000  psi 
and  arc  plotted  as  a  single  curve  on  Fig.  29.  The  theo¬ 
retical  curves  in  Fig.  29  arc  plotted  assuming  a  value  of 
N'/>  =  2  which  corresponds  to  a  value  of  £>j  S  20°.  The 
results  of  the  direct  shear  tests  shown  in  Fig.  7  show 
that  the  value  of  <>j  decreases  with  increasing  normal 


Figure  28.  Stress-sti  on  curve  used  for  plexiglass  liners. 
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pressure  on  the  joints.  The  direct  shear  tests  indicated 
that  between  a  normal  pressure  of  SO  psi  and  400  pst, 
the  value  of  decreased  from  33“  to  29“.  The  r 
pressures  of  interest  ranged  up  to  S00  psi  and  the  ma¬ 
terial  around  the  opening  would  also  be  subjected  to  a 
stress  concentration  due  to  the  presence  of  the  tunnel. 
Thus  it  was  decided  to  compare  the  experimental  re¬ 
sults  of  the  nit  del  tests  to  theoretical  curves  calculated 
using  a  value  of  </>j  s  20°.  In  big.  29,  the  experimental 
curve  starts  out  well  above  the  theoretical  curves  indi¬ 
cating  that  the  effective  value  of  0j  at  low  pressure  was 
considerably  higher  than  20°  as  should  be  expected.  At 
a  model  pressure  of  400  psi,  the  experimental  curve 
crosses  the  theoretical  curve  for  N</>  =  2.  o„  =  200  psi. 
With  increasing  model  pressure,  the  experimental  curve 
crosses  over  the  theoretical  curve  for  N0=  2.  ou  =  100 
psi.  and  then  becomes  parallel  to  the  theoretical  curve 
for  N  0  =  2,  au  =  0.  Thus  the  ■phenomenon  occurring  in 
the  models  could  be  considered  to  be  a  decreasing  ef¬ 
fective  unconfined  strength  of  the  jointed  mass  with 
increasing  model  pressure  at  constant  N</>.  As  a  matter 
of  fact,  the  value  of  N  also  decrea  es  with  increasing 
model  pressure  and  a  more  detailed  analysis  would  take 
this  effect  into  account  quantitatively.  However  the 
model  material  properties  used  arc  estimates  and  so¬ 
phisticated  refinements  of  the  analysis  arc  not  justified. 

Fig.  30  is  a  plot  similar  to  Fig.  29  with  the  experi¬ 
mental  results  of  two  model  tests  (JB  *7  and  JB  —1 0) 
shown  with  theoretical  curves  calculated  from  assumed 
values  of  the  model  properties.  JB  #7  had  two  sets  of 
joints  spaced  at  I  in.  and  oriented  at  45°  to  the  princi¬ 
pal  loading  directions  and  had  .,  4-in.  diameter  tunnel 
lined  with  a  1/8-in.  thick  plexiglass  liner.  JB  #10  had 
two  sets  of  joints  spaced  at  1  in.  and  oriented  parallel 
to  the  principal  loading  directions  with  an  <S-in.  diame¬ 
ter  tunnel  lined  with  a  I  /4-in.  thick  plexiglass  liner. 
These  two  liners  have  the  same  calculated  circumfer¬ 
ential  stiffness  (Et/R  =  25,000  psi)  assuming  the  plexi¬ 
glass  was  the  same  in  each  liner.  The  average  measured 
model  stiffness  of  JB  #7  (Em  =  64,300  psi)  was  greater 
than  that  of  JB  #10  (Em  =  61,000  psi) 

The  experimental  results  of  JB  #7  and  JB  #10  plot 
very  nearly  on  the  same  line  up  to  a  model  pressure  of 
about  500  psi  above  which  they  diverge  rather  sharply. 
The  experimental  curve  for  JB  #10  starts  out  above  all 
the  theoretical  curves  and  then  at  a  model  pressure  of 
about  300  psi  becomes  coincident  with  the  experi¬ 
mental  curve  for  N0  =  2,  «„  =  300  psi  and  follows  that 
curve  closely  up  to  a  model  pressure  of  950  psi.  There, 
it  begins  to  rrovc  toward  curves  lor  lower  values  of  rru. 


lhc  experimental  curve  for  JB  #7  also  starts  out  above 
.ill  the  theoretical  curves,  but  it  then  crosses  over  all 
the  experimental  curves  between  model  pressures  of 
400  psi  and  550  psi. 


5  CONCLUSIONS 

The  results  of  five  model  tests  arc  reported  and 
analv/ed  in  tins  repoit.  T  he  tunnel  liners  used  in  these 
models  were  chosen  so  that  the  ratio  of  the  circumfer¬ 
ential  tunnel  liner  stiffness  to  rock  mass  stiffness 
(Et/R)/iim  was  similar  to  that  for  concrete  liners  in 
rock.  Tin-  plexiglass  liners  used  in  this  study  had  the 
further  advantage  that  the  ratio  of  the  circumferential 
stiffness  iF.i/R)  to  the  flexural  stiffness  ( Fl/R* )  was 
nearly  the  same  as  tor  reinforced  concrete  tunnel  liners. 
In  the  model  tests  reported  here,  none  of  the  plexiglass 
liners  buckled  and  twoof  them  failed  in  circumferential 
compression  :.s  a  real  reinforced  concrete  liner  wor  ld. 
Thus  t  ie  use  ot  plexiglass  liners  prevented  a  buckling 
mode  of  failure  which  had  previously  been  experienced 
with  aluminum  liners.1* 

The  analysis  of  the  data  reported  herein  shows 
that  elastic  theory  is  not  sufficient  to  predict  the  be¬ 
havior  of  models  loaded  to  high  model  pressures.  In 
these  tests,  elastic  theory  was  sufficiently  accurate  for 
calculating  the  diametrical  strains  of  the  lined  tunnels 
up  to  diametrical  strains  of  about  1.2%.  For  diametri¬ 
cal  strains  greater  than  1.2%,  calculations  from  clastic 
theory  underestimate  the  diametrical  .drains  of  a  liner. 
An  elusto-plastic  theory  developed  by  Hcndron  and 
Aiyer  was  used  to  analyze  the  model  behavior  in  the 
plastic  range.  Using  experimentally  measured  strength 
and  stiffness  parameters  for  the  model,  the  Hendron- 
Aiyer  analysis  could  be  used  to  estimate  the  diametri¬ 
cal  strains  of  the  tunnel  liner  in  the  plastic  range  of 
behavior. 

The  detailed  analysis  of  the  experimental  data 
presented  in  the  previous  chapter  shows  that  the  pro¬ 
cedure  described  below  can  be  used  for  prediction  of 
tunnel  liner  deformations  (AD/D.  diametrical  strains) 
for  a  tunnel  liner  in  a  jointed  rock  mass  subjected  to 
hydrostatic  compression.  The  Young’s  modulus  of  the 
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jointed  rock  mass,  Em,  as  observed  in  these  tests  was 
about  1/1  Oth  of  the  Young’s  modulus.  E,  of  the  intact 
model  material  which  composes  the  joint  blocks.  If  the 
Young's  modulus  of  the  jointed  mass,  Em.  is  used  in  an 
elastic  analysis  of  the  liner  and  jointed  mass,  it  was 
found  that  the  deformation  of  the  liner  would  be  pre¬ 
dicted  within  about  ±  10%  for  diametrical  strains 
(AD/Dj  less  than  1.2%.  For  strains  above  1.2%.  the 
elas to- plastic  analysis  given  by  Hendron  and  Aiycr 
could  be  used  to  predict  the  liner  deformations  if  the 
jointed  medium  was  assumed  to  have  the  following 
properties: 

( 1)  The  Young’s  modulus  of  the  jointed  medium. 
Em,  in  the  elastic  range  should  be  taken  as 
about  I /I Oth  the  Young's  modulus  of  the  in¬ 
tact  model  material. 

(2)  The  angle  of  friction  of  the  jointed  medium. 


0m.  should  be  taken  as  the  angle  of  shearing 
resistance  along  the  joints,  <?j. 

(3)  The  effective  unconfined  strength  of  the  joint¬ 
ed  medium  surrounding  the  tunnel,  ou,  should 
be  taken  as  a  fraction  of  the  unconfincd 
strength  of  the  intact  rock  material  as  shown 
in  Fig.  27,  the  fraction  decreasing  with  an  in¬ 
crease  in  the  ratio  of  tunnel  diameter  to  joint 
spacing  as  shown  in  Fig.  27. 

If  the  above  procedure  is  followed,  the  diameter 
changes  of  the  structural  liners  can  be  predicted  with¬ 
in  ±  1 5%  for  the  test  results  piesentcd  in  this  report.  It 
has  been  the  writers’  experience  that  the  same  pro¬ 
cedure  can  be  applied  to  concrete  tunnel  liners  in  rock. 
But.  if  this  procedure  is  used  for  concrete  liners  in 
rock,  clastic  analyses  should  not  be  used  beyond  a  dia¬ 
metrical  liner  strain  of  about  0.3%c 
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APPENDIX  A:  SIMILITUDE  CONSIDERATIONS 


Fundamental  Considerations.  In  most  physical  phe¬ 
nomena  considered  in  civil  engineering  and  geology  ir  is 
conventionally  assumed  that  a  “cause  and  effect"  re¬ 
lationship  exists  between  the  various  independent  and 
dependent  variables  which  influence  and  describe  a 
phenomenon.  This  relationship  is  assumed  to  be  ex¬ 
pressed  by  some  function 

f(X|,  xj,  x,.. ...x . xn)=  0  | Eq  1 1 

where  Xj  are  pertinent  independent  and  dependent  vari¬ 
ables.  The  function  f<Xl)  can  usually  be  expressed  ex¬ 
plicitly  for  only  the  simpler  phenomena.  It  may  be 
determined  from  either  theoretical  considerations  or 
empirical  studies.  Because  the  most  basic  physical  laws 
(such  as  Newton’s  Laws)  are  dimensionally  homogene¬ 
ous;  that  is,  their  form  does  not  depend  upon  the  units 
of  measurement,  it  can  be  hypothesized  that  the  more 
complex  function  f<x.(  is  also  dimensionally  homo¬ 
geneous,  even  though  it  is  not  explicitly  known. 

The  theory  of  dimensional  analysis,  founded  in  the 
mathematical  theories  of  algebra,  is  summarized  in 
Buckingham’s  theorem,29  which  essentially  states  that 
from  the  dimensionally  homogeneous  function  or 
equation  describing  a  phenomenon,  it  is  possible  to 
develop  a  relationship  in  which  the  variables  appear  in 
a  set  of  dimensionless  products.  (For  a  more  complete 
discussion  see  texts  such  as  Murphy  or  Langhaar.)  In 
practice,  dimensional  analysis  allows  us  to  determine 
these  dimensionless  products,  given  the  pertinent  vari¬ 
ables,  even  though  we  do  not  know  the  form  of  the 
function  f(Xj)  which  describes  the  phenomenon.  Thus, 
from  Equation  1,  which  is  the  basic  function  relating 
the  pertinent  variables  x,  in  a  description  of  the  phe¬ 
nomenon,  we  arrive  by  a  dimensional  analysis  at 

F  (jr,  ,7T;  ,7r3 . 7Tj . rrm )  =  0  |Eq2| 

as  a  description  of  the  phenomenon.  Each  term  7Tj. 
often  called  a  Buckingham  pi  term,  is  a  dimensionless 
product  of  some  number  of  the  original  Xj  variables. 
Generally,  the  number  m  of  independent  pi  terms  is 
related  to  the  number  n  of  the  xj  variables  and  to  the 
number  r  of  fundamental  dimensions  (such  as  mass, 
length,  time,  temperature)  which  arc  involved  in  the  x, 
variables  by 

m  =  n  r  |Eq  3] 


”11.  L.  Langhaar.  Dimensional  Analysis,  p  18 


There  are  several  advantages  which  may  be  gained 
from  the  dimensional  analysis.  First,  the  relationship 
between  the  sj  variables  in  the  Wj  terms  often  gives 
valuable  insight  into  the  phenomenon  being  con¬ 
sidered.  Second,  the  phenomenon  is  described  in  terms 
of  a  fewer  number  of  variables,  only  m  pi  terms  instead 
of  the  n  original  variables.  This  reduction  of  variables  is 
often  important  when  studying  the  phenomenon  ex¬ 
perimentally  since  it  usually  reduces  the  number  of 
experiments  which  must  be  conducted.  A  third  advan¬ 
tage  is  that  the  dimensional  analysis  provides  a  theo¬ 
retical  basts  for  model  studies,  by  which  it  may  be 
possible  to  reduce  even  further  the  cost  and  time  in¬ 
volved  in  studying  the  phenomenon  being  considered. 

The  function  given  in  Equation  2  is  dimensionally 
homogeneous  and  completely  general.  If  the  pi  terms 
are  independent  and  contain  all  of  the  pertinent  vari¬ 
ables  Xj  which  influence  and  describe  the  phenomenon, 
then  the  function  Fj^.)  completely  describes  the  phe¬ 
nomenon.  regardless  of  the  scale  of  units  with  which 
the  quantities  x;  are  measured,  and  regardless  of  the 
absolute  magnitude  of  the  x,  quantities.  This  means 
that  if  we  wish  to  utilize  models  to  study  the  behavior 
of  a  prototype,  we  can  he  assured  that  the  behavior  of 
the  model  duplicates  the  behavior  of  the  prototype  in 
a!i  respects  if  each  of  the  pi  tcims  for  the  model  is 
equal  to  the  equivalent  pi  term  for  the  promt-  pc:  that 
is.  if 

(Th)  model  ~('u)  prototype  |Eq4] 

If  such  a  condition  exists,  all  requirements  of  simili¬ 
tude  have  been  satisfied,  the  model  is  said  to  be 
“completely  similar”  to  the  prototype,  and  the  phe¬ 
nomenon  in  the  model  is  an  exact  replica  of  the  phe¬ 
nomenon  in  the  prototype. 

Several  major  problems  develop  in  the  practical 
application  of  dimensional  analysis  to  modeling.  One  is 
that  it  generally  is  technically  impossible  to  insure  that 
all  of  the  Xj  variables  which  influence  the  prototype  are 
con.  Jered  in  the  dimensional  analysis  and  arc  accu¬ 
rately  reproduced  in  the  model  by  equating  the  t q 
terms  of  the  model  and  prototype.  For  satisfactory 
modeling  it  is  necessary  that  the  phenomenon  be 
understood  well  enough  to  know  what  variables  (xj) 
are  most  significant  and  which  pi  terms  (rq)  must  be 
duplicated  most  rigorously  so  that  the  model  gives  the 
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nuiM  accurate  simulation  of  the  prototype  which  it  is 
both  possible  and  practical  to  attain. 

In  this  approximation  of  the  prototype  by  the 
model  a  second  major  problem  arises  from  what  are 
known  as  “scale  effects.”  For  example,  as  physical  size 
o!  the  model  varies,  the  relative  importance  of  differ¬ 
ent  forces  may  vary  also.  Body  forces  such  a-  weight 
due  to  gravitational  attraction  vary  as  the  mass  of  the 
body,  hence  as  the  third  power  of  linear  dimensions, 
while  surface  forces  such  as  pressures  vary  as  area, 
hence  as  the  second  power  of  linear  dimensions.  Thus, 
as  the  physical  size  of  a  model  is  reduced,  the  influence 
of  body  forces  decreases  more  rapidly  than  that  of 
surface  forces.  It  is  possible  that  the  hehavroi  of  the 
prototype  may  be  strongly  influenced  In  body  forces, 
but  when  a  model  is  made  at  a  reduced  size,  its  be¬ 
havior  may  be  strongly  influenced  by  surface  forces 
which  arc  of  lesser  importance  in  die  prototype,  fate¬ 
ful  consideration  must  be  given  n>  problems  such  as  the 
selection  of  pertinent  variables  and  the  possibi’  !y  ot 
scale  effects  in  developing  a  model  testing  program. 

Selection  of  Significant  Variables.  The  dimensions  ot 
mass  M,  length  L.  and  time  (  are  probably  the  most 
commonly  used  dimensions  by  •..»•  '  physical  phe¬ 
nomena  are  described.  An  c<|iia!ly  valid  set  of  basic 
dimensions  is  force  F,  length  L,  and  time  T.  For  a 
static  system  such  as  is  being  considered,  only  force  F 
and  length  L  are  involved.  This  is  the  set  which  will  be 
used  in  the  following  analysis. 

The  significant  variables  associated  with  the  be¬ 
havior  of  an  underground  opening  which  will  be  con¬ 
sidered  in  this  study  are  given  in  Table  Al.  The  reasons 
for  the  selection  of  these  variables  and  the  exclusion  of 
others  is  discussed  below.  A  number  1  in  the  dimen¬ 
sions  column  means  that  the  variable  is  dimensionless, 
a  pure  member. 

Frcc-fieU  stresses.  The  prototype  chosen  for 
study  is  a  segment  of  a  long,  straight,  horizontal  tunnel 
buried  ut  derground  at  a  depth  several  times  greater 
than  the  tunnel  diameter.  The  most  significant  forces 
influencing  the  tunnel  behavior  are  assumed  to  he 
those  due  to  the  frcc-field  stresses  which  would  exist  at 
the  location  of  the  tunnel  if  it  were  not  present.  The 
two  free-flcld  stresses  considered.  <\  and  otl.  are  the 
vertical  and  horizontal  stresses  in  a  plane  perpendicular 
to  the  tunnel  axis  (Fig.  Al),  and  are  assumed  to  be 
principal  stresses  (which  they  would  be  in  an  elastic 
half  space  with  a  horizontal  surface).  The  magnitude  of 
the  third  principal  stress,  uy.  the  horizontal  stress  paral¬ 
lel  to  the  tunnel  axis,  will  he  considered  later. 


I  he  veitical  fiee-fieki  stress  uv  could  be  due  to  the 
weight  of  the  overlying  material,  and  at  any  depth  z 
below  the  surface  it  would  be  given  by 

jF.q  5 1 

where  7  is  the  average  unit  weight  of  the  overburden. 
Following  the  arguments  of  Ter/.aghi  and  Richart.  ,n 
and  Deere,"  the  hoii/onta!  free-flcld  stress  <>1,  is  as¬ 
sumed  to  be  given  by 

oh  =  N'ov  =  N  7  z  |  F.q  6 1 

file  coefficient  N  (dating  uv  and  oj,  is  intimately  re¬ 
lated  to  the  present  geologic  environment  and  the  pre¬ 
cious  geologic  history  of  the  site  and  rnay  vary  over  a 
wide  range  of  values. 

The  tree-held  stresses  could  also  be  due  to  the 
loads  imposed  by  a  nuclear  detonation  over  a  protec¬ 
tive  Ntou  lute  in  rock.  In  this  case  the  major  principal 
sfess  would  be  the  dynamic  radial  stress.  (iri|.  of  the 
c’iicct  induced  ground  shock  emanating  from  the  era- 
let.  The  minor  ptmcqra!  stress  would  be  the  dynamic 
tan centia!  stress.  given  by 

•W  =  «rj(yjy)  I1 -A  7 1 

where  v  is  Poisson's  ratio  of  the  rock  mass.  The  ratio 
between  the  major  and  minor  principal  stresses.  N.  for 
this  case  could  range  from  about  I  /4  to  1.0. 

!<■  -arise  the  value  ol  N  for  both  static  and  dy¬ 
namic  si.  Acids  varies  over  a  large  range,  model  tests 
to  simulate  either  of  these  problems  should  be  con¬ 
ducted  at  vat  ions  values  of  N  rather  than  at  a  specific 
value  of  N  lints  lor  purposes  ol  dimensional  analysis 
of  the  static  problem  considered  here  both  <jv  and  «q 
will  be  considered  as  independent  variables. 

Miiidliir'1 1  has  developed  rigorous  solutions  f  ir 
stress  distributions  around  a  circular  tunnel  in  an  elas¬ 
tic  half  space  under  gravity  loading  with  N  values  r»f  0. 
pits  ami  i.  His  results  (see  Pattck35  or  Camlh.  and 

50  ts.  Icr/:i"hi  and  I  i  Richart  Jr.,  "Stresses  ill  Rod.  \bollt 
<  avitjc.,”(.'e«>/e(7u;j^»c,  Vo!  i  1 19.52)  pp  57-9(1. 

'MM1  Deere  Diseuscion  of  •’failure  of  llomop.'neo.is  Rock 
('ruler  Dynamic  Compressive  la  rad  in!'.”  try  (i.Ii.  Clark  and 
R  D  (  audit-.  Sene  "I  Stress  in  the  Earth's  Crust.  W.lt.  Judd, 
ed.  !l  Isevier.  I9f,4|  p;r  321-323;  ami  "Geologii  Considcr- 
.itionv.*’  Rink  Metli.  in  Eng.  Praet..  Stapj;  and  Zienkiewic •/. 
ed.  Irtlin  Wilcv  .V  Sons.  19frK). 

oR.D.  Mindlin.  "Stress  Distribution  Around  A  Tunnel,”  Prtte. 
I.SY  .’  r  April  DM'*)  pp  M‘l  n42 

rr|  ,.\  patiek .  Stresses  -Mu nit  Stine  Openings  in  a  Homogene- 
•  «»-.  /.’or  k  liinlv  I  New  York.  1451  t 
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Table  A I 

Significant  Variables 


Variable 

Dimensions 

ITce-field  Slfcs-.es,  assumed  to  be  principal  stresses 

°v 

the  vertical  tree-held  stress 

1  I.2 

"to 

the  horizontal  tree-field  stress 

11.-’ 

Intact  rock  properties 

conn,. 

cohesion  or  uncontined  compressive  stiength.  -t-e  either  one 

II.-1 

0. 

angle  of  internal  friction 

1 

tensile  strength 

II/2 

c. 

modulus  of  elasticity 

IV2 

U, 

Poisson’s  ratio 

1 

Rock  mass  properties 

s. 

spacing  of  |oints 

1. 

Oj. 

angle  of  frictional  resistance  along  joint  surfaces 

1 

n% 

orientation  of  joint  planes  with  respect  to 
principal  planes 

1 

Geometry 

of  the  opening 

d  or  a. 

internal  diameter  or  radius  of  unlined  opening, 
and  1)  or  R.  diameter  or  radius  of  the  lining 

1. 

Support  properties 

I't  l/RJ- 

flexural  stiffness  of  lining,  where  l  g  is  the 
modulus  of  the  liner  material  anil  1  is  the  moment 
of  inertia  of  the  liner  wall  per  unii  length 

IT/2 

I'iyt/R, 

compression  stiffness  of  lining,  wlvrre  t  is 
the  lining  thickness 

II/2 

Response  of  system 

u. 

'adial  deformations  of  the  opening  wall 

I. 

f. 

strain  within  the  rock  mass  around  the  opening 

1 

o. 

stress  within  the  rock  mass  around  the  opening 

II/2 

r. 

thrust  in  the  lining,  force  per  unii  length 

II/* 

M, 

moment  in  the  lining,  moment  per  unit  length 

I- 

Clark34  for  more  detail)  show  that  if  the  tunnel  is  at  a 
depth  of  three  tunnel  diameters  or  more,  the  stress 
distribution  around  the  tunnel  can  be  approximated 
very  closely  by  the  distribution  of  stresses  about  a  cir¬ 
cular  hole  in  a  biaxially  loaded  clastic  flat  plate,  as 
determined  by  the  Kirsch  equations  (see.  for  example. 
Timoshenko  and  Goodior.3 s  or  Obert  and  Duval3'1  i.  In 


J4R.D.  Caudle  and  Ci.li.  ('lark.  “Stress  Around  Mine  Openings 
in  Some  Simple  Geologic  Structures.”  University  nf  Illinois 
Hug  Exp.  Sta.,  Bull.  No.  4M)  (1955). 

s»S.  Timoshenko  and  J.N.  Goodior,  Theory  of  Elasticity 

(McGraw-Hill,  1951). 

1 6 1..  Obert  and  W.  I.  Duvall,  Rock  Mechanics  am I  the  Design 
of  Structures  in  Rock.  (John  Wiley  &  Sons,  I9(i7). 


both  cases  the  area  in  which  stresses  are  significantly 
influenced  by  the  tunnel  has  a  width  of  about  four 
times  the  diameter  of  the  tunnel. 

These  observations  allow  a  considerable  simpli¬ 
fication  in  modeling  the  stress  field  about  the  tunnel. 
The  observation  that  the  area  influenced  by  the  tunnel 
is  about  four  tunnel  diameters  wide  means  that  the 
stress  distribution  on  a  square  whose  sides  are  four 
tunnel  diameters  long,  concentric  with  the  tunnel,  is 
very  closely  approximated  by  that  illustrated  in  Fig. 
A I  where  ov  and  are  the  frcc-ficld  stresses  at  the 
location  of  the  tunnel  axis,  and  &nv  and  ^0|,  are  the 
changes  in  the  frce-field  stresses  between  the  top  and 
the  bottom  of  the  /.one  ilue  to  the  weight  of  the  ma- 
terial  within  it. 
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Figure  Al.  Stress  distribution  some  distance  from  tun¬ 
nel. 


Figure  A2.  Approximate  stress  distribution  some  dis¬ 
tance  from  tunnel. 


The  observation  that  Mindlin’s  solutions  for  tun¬ 
nels  buried  at  more  than  three  tunnel  diameters  are 
closely  approximated  by  Kirsch’s  solution  means  that 
the  stress  distribution  of  Fig.  A I  can  be  approximated 
by  that  of  Fig.  A2  with  small  error.  This  means  that  for 
the  clastic  case  the  stress  distribution  around  a  tunnel 
underground  at  a  depth  of  more  than  three  tunnel 
diameters  is  more  strongly  influenced  by  the  average 


free- field  stresses  at  the  tunnel  location  than  by  the 
gravity  forces  on  the  material  within  the  zone  influ¬ 
enced  by  the  tunnel.  When  modeling  the  static  behav¬ 
ior  of  a  tunnel  in  this  situation  only  the  average  free* 
field  stresses  as  illustrated  by  Fig.  A2  need  be  modeled 
and  the  body  forces  can  be  neglected.  On  the  basis  of 
the  preceding  arguments  the  body  forces  and  associ¬ 
ated  variables  such  as  the  density  of  the  material  arc 
not  considered  in  the  list  of  significant  variables  for  the 
present  study.  Thus,  forces  due  to  the  weight  of  loos¬ 
ened  rock  around  the  opening  are  not  being  modeled. 

Intact  rock  properties.  A  discussion  of  the  failure 
mechanism  of  rock  materials  is  beyond  the  scope  of 
this  report.  The  reader  is  referred  to  papers  such  as 
Jaeger,17  Hieniawski,38  Hendron,39  and  Hock40  as 
examples  and  discussions  of  recent  work  in  this  area. 
Regardless  of  the  actual  failure  mechanism  in  rock 
materials,  it  is  generally  observed  that  some  form  of 
Mohr  envelope  can  lx  fitted  to  observed  experimental 
data  and  used  to  predict  rock  strength.  In  order  to 
simplify  the  dimensional  analysis  and  subsequent  dis¬ 
cussions,  it  will  be  assumed  that  the  general  curvilinear 
Mohr  envelope  can  he  approximated  by  a  straight  line 
in  the  compression  pressure  range  of  interest,  so  that 
the  general  Mohr  failure  criteria 

r  -  F(0)  I  Be*  8) 

can  he  replaced  by  the  more  specialized  Cottlomb- 
Navier  failure  criteria 

r  =  c  +  a  tan  <>  |  Eq  0] 

(For  a  more  detailed  discussion  of  failure  criteria  sec 
Nudai.41  Seely  and  Smith,43  Jaeger,43  and  Oberl  and 
Duval.) 

The  two-independent  Coulomb-Navicr  strength 


5  ‘  J.C.  Jaeger,  “Brittle  Fracture  of  Rocks,"  Sill  Svmp.  on  Rock 
Merli..  (A IMF,  1967)  pp  3-131. 

5  “Z.T.  Btenunvski,  “Mechanism  of  Brittle  Fracture  or  Rock.” 
Ini.  Jotirn,  Rock  Mfch.  and  Mining  Sci.  (October  1967)  pp 
3954  30. 

J,A.J.  Hcmlron.  Jr.,  “Mechanical  Properties  of  Rock,"  Rock 
Merit,  in  ling.  I’m  cl..  Staae  &  Zienkicwiez.  ed.  (John  Wilev 
&  Sons.  1968). 

I  Hock.  ''Brittle  Failure  of  Rock.” 

* 1  A.  Njil.it.  Thcor,’  of  How  ami  Fracture  of  Solids,  Vols  I  ami 

II  (McGraw-Hill,  1950). 

47  F.B.  Seelv  ami  J.O.  Smith,  Advanced  Mechanics  of  Materials 
(John  Wiley  A  Sons.  I  '>52). 

41  J.C.  Jaeger.  I  laslicilv,  Fracture  and  Flow  (Mctliucn  and  Co. 
I  td..  I *>*.21. 
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parameters  are  the  cohesion  c  and  the  angle  of  internal 
triction  0.  An  alternate  and  equally  valid  pair  are  the 
unconfined  compressive  strength  qu,  and  the  angle  of 
internal  triction  0.  Either  r  or  q„,  in  conjunction  with 
0.  arc  necessary  and  sufficient  to  define  the  failure 
state  and  describe  the  failure  envelope  in  the  compres¬ 
sion  range.  In  addition,  the  tensile  strength  ot  must  be 
defined,  since  Equation  ‘)  is  not  valid  in  the  tensile 
stress  range.  The  complete  failure  envelope,  then,  is  as 
shown  in  Fig.  A3.  A  rigorous  consideration  of  the  actu¬ 
al  curvilinear  Mohr  envelope  does  not  change  the  basic- 
conclusions  of  the  dimensional  analysis. 

The  elastic  constants,  the  modulus  of  elasticity  E 
and  Poisson’s  ratio  u,  relate  stress  and  strain  assuming 
the  intact  rock  exhibits  a  quasi-elastic  behavior  at  low 
and  intermediate  stress  levels.  In  general  this  quasi- 
elastic  range  is  followed  by  a  range  in  which  inelastic 
strains  occur,  and  then  by  some  form  of  failure  such  as 
fracture  or  plastic  deformation  (Fig.  A4).  No  variables 
are  included  to  describe  the  inelastic  and  plastic  regions 
for  two  reasons:  ( 1 )  because  of  the  wide  range  of  be¬ 
havior  exhibited  by  different  rocks,  and  (21  because  of 
the  scarcity  of  real  numerical  description  and  data  tor 
this  portion  of  the  stress-strain  curve.  It  is  recognized 
'•“.at  the  other  types  of  behavior  exist.  For  example,  a 
cot.cavc  upward  stress-strain  curve  is  commonly  ob¬ 
served  at  low  stress  levels  for  very  porous  rocks,  for 
highly  weathered  i  ticks,  and  for  thinly  bedded  ot  foli¬ 
ated  rocks  compressed  perpendicular  to  the  bedding  or 
foliation.  The  initial  quasi-elastic  behavior  is  possibly 
more  common,  however,  and  is  much  more  simple  to 
consider  and  model. 

Time-dependent  behavior  such  as  creep  or  viscous 
deformation  is  not  considered  here.  These  properties  of 
real  rock  are  so  poorly  known  and  understood  that  any 
attempt  to  consider  them  in  modeling  the  behavior  of 
underground  openings  must  be  considered  a  very  ques- 
ionable  practice  for  the  “present  state  of  the  art."  The 
ore  exception  to  this  statement  would  be  in  the  case  of 
underground  openings  in  the  evuporites:  rock  salt,  pot¬ 
ash,  and  possibly  gypsum  and  anhydrite.  For  these 
rocks  time-dependent  behavior  is  so  pronounced  (hat  it 
dominates  the  behavior.  Indeed,  because  it  is  so  pro¬ 
nounced  it  can  he  and  has  been  studied  enough  so  that 
intelligent  attempts  to  model  the  time-dependent  be¬ 
havior  of  such  rocks  can  be  made  (see  for  example 
Thompson  and  Rippergcr44 ). 


44  li.  Thompson  and  K.A.  Rippergcr  "An  Kxpcrinientcd  Tech¬ 
nique  for  the  Investigation  of  the  flow  of  Halite  and  Syl- 
vanite,"  Sixth  Symposium  on  Rock  Mechanics  (University 
of  Missouri,  1964)  pp  467488. 


Rock  Mass  1‘ropcrties.  The  most  important  rock 
mass  properties  of  a  rock  mass  surrounding  a  tunnel  are 
the  spacing  and  oiicntution  of  the  discontinuities  and 
the  strength  along  the  discontinuities.  Thus  for  pur¬ 
poses  of  dimensional  analysis,  the  spacing  of  the  joints 
or  discontinuities,  s.  should  he  considered  as  an  inde¬ 
pendent  variable.  The  orientation  of  each  joint  set  in  a 
particular  problem  is  also  very  important  and  must  be 
considered  as  an  independent  variable.  It  is  suggested 
herein  to  specify  the  orientation  of  the  joint  planes,  0 
with  respect  to  the  maximum  free-field  principal  stress 
direction  since  the  orientation  will  govern  the  ratio  of 
normal  and  shear  stresses  on  the  joint  planes. 

The  strength  properties  of  the  joint  surfaces  must 
also  he  considered.  In  the  model  tests  conducted  on 
this  study  the  models  were  composed  of  planar  joint 
surfaces  free  of  irregularities.  The  shear  strength  along 
such  a  plane  surface  can  be  expressed  by 

m  -  tan  0j  |F.q  1 01 

where  o„  is  the  normal  stress  on  the  plane  surface  and 
is  a  property  of  the  joint  surfaces  which  should  be 
included  in  the  dimensional  analysis  to  represent  the 
angle  of  shearing  icsisiancc  along  the  planar  joint  sur¬ 
faces. 

Opening  geometry.  For  this  study,  the  opening  is 
assumed  to  be  of  circular  shape  with  internal  diameter 
d  (radius  a).  It  is  assumed  that  any  lining  of  the  open¬ 
ing  is  of  diameter  D  (radius  R)  which  is  the  same  mag¬ 
nitude  as  the  opening  diameter  d.  and  that  the  variables 
d.  a.  D.  and  R  may  be  used  interchangeably  for  the 
purposes  of  dimensional  analysis.  However,  only  one  of 
these  variables  may  be  considered  as  an  independent 
variable  for  the  purposes  of  dimensional  analysis,  since 
they  ail  •‘•.•fine  the  same  thing,  the  size  of  the  under¬ 
ground  opening. 

Lining  properties.  “Unsupported”  cylinders  and 
rings:  general  expressions  for  the  clastic  deflection  of 
rings  and  cylinders  due  to  loads  which  do  not  vary  with 
position  along  the  axis  arc  developed  in  numerous  texts 
on  strength  of  materials,  the  theory  of  elasticity,  and 
structural  analysis  (stn.it  as  Seclv  and  Smith  or  Timo¬ 
shenko  and  Gere45  ).  It  is  seen  from  these  expressions 
that  the  deflection  of  a  point  on  a  long  cylinder  due 
to  bonding  deformation  caused  by  external  loads  P  is 
given  by 

5b=C,  PR1  (lV)/E«l  | Eq  11] 


4,S.  Timoshenko  amt  J.  (icre.  Theory  of  plastic  Stability 
(Me<ira\v-Hi)l.  1961). 
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Figure  A3.  Failure  envelope  considered  typical  for  rock. 
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and  the  deflection  5C  due  to  circumferential  compres¬ 
sion  deformation  is  given  by 

&c=Cj  PR(1V)/Evt  (Eq  12J 

Tlie  quantities  C(  and  Cj  are  constants  whose  values 
depend  upon  the  i  omtguration  of  the  applied  loads 
and  upon  the  location  of  the  point  in  question.  Lane46 
and  Dorris47  give  tables  of  values  of  C,  for  different 
loading  configurations.  The  flexural  and  compression 
stiffnesses  of  a  long  cylinder  are  the  quantities  E  I/R3 
and  Evt/R(  1  -t>v * )  respectively.  The  term  (I 
•i>v2 )  arises  because  of  the  restraint  offered  by  the  axial 
stress  parallel  to  the  cylinder  axis  in  plane  strain  prob¬ 
lems.  In  the  case  of  a  ring  subjected  to  plane  stress,  this 
term  vanishes  and  the  bending  stiffness  becomes 
Ejl/R3. 

Ruckling  of  unsupported  rings  and  cylinders  under 
uniform  and  non-uniform  external  loads  and  pressures 
is  treated  in  a  number  of  texts  and  papers  such  as  Seely 
anil  Smith,  Timoshenko  and  Gere.  Boresi  4  R  Bodner 4  9 
and  Anderson  and  Boresi.30  The  conclusion  reached  in 
these  investigations  is  that  for  a  long,  thin-walled  cylin¬ 
der.  the  critical  external  buckling  pressure  pcr  is  given 
bv  an  equation  of  the  form 

p,r  =  C3  Evl/R3(lV)  |EqU| 
where  C3  is  a  function  of  the  loading  configuration  and 
the  buckling  mode.  Note  that  the  significant  cylinder 
parameter  is  the  bending  stiffness  E  l/(l*Uv2)R3,  which 
is  reasonable  since  buckling  is  a  bending  phenomenon. 
As  before,  the  term (1 -up2)  is  due  to  restraint  from  the 
third  dimension. 

Cylinders  and  rings  on  clastic  foundations:  Helen- 
yi51  shows  that  for  a  ring  on  elastic  foundation  (i.e.. 


**  K.S.  Lane,  “Garrison  Darn  Test"  Tunnel,  Evaluation  of  Re¬ 
sults,"  Trans.  ASCF.  Vol  125.  Pt  I.  Paper  3022  (I960)  pp 
268-306. 

4  ’  A.I-.  Dorris.  Response  of  Horizontally  Oriented  Buried 
Cylinders  to  Static  and  Dynamic  Loading,  Technical  Report 
No.  AI'WL-TR-65-l  16  (Air  force  Weapons  Laboratory, 
1966). 

44  A. P.  Uorcsi,  “A  Refinement  of  the  Theory  of  Buckling  of 
Rings  Under  Uniform  Pressure.”  Journal  of  Applied  Mechan¬ 
ics.  Vol  22  (ASML,  1955). 

4#S.R.  Bodner,  “On  the  Conservative  ness  of  Various  Distri¬ 
buted  force  Systems.”  Journal  of  the  Aeronautical  Sciences. 
Vol  25.  No.  2  (1958). 

50  R.ll.  Anderson  and  A.P.  Boresi.  “Equilibrium  and  Stability 
of  Rings  Under  Nonunilormly  Distributed  Loads."  Proc.  of 
the  Fourth  (J.S.  Mat.  Cong,  of  Applied  Merit.,  Vol  I  (ASMf, 
1962). 

’ 1  M.  Ilelenyi,  Beams  on  Clastic  Foundations  (University  of 
Michigan  1946). 


with  radial  spring  supports)  which  is  subjected  to  any 
generalized  tadial  loading  P,  the  radial  displacement  5 
of  any  point  of  the  ring  is  of  the  form 

6  =  F<P/t:pl/R\  KR/Kpl/R3 ,  o)  [Eq  141 

where  K  is  the  modulus  of  subgrade  reaction  (the 
“soring  constant”  of  the  foundation),  and  0  is  a  meas¬ 
ure  of  position  on  the  ring.  In  the  case  cf  a  long  cylin¬ 
der  instead  of  a  ring,  the  term  Eel/(t-ue) is  substituted 
in  place  of  EpI.Note  that  this  is  (l-uB),  not  (l-uB2)as 
previously,  the  reason  being  that  the  elastic  foundation 
reaction  is  not  continuous  but  consists  of  discrete 
springs. 

Cheney51  and  Luscher53  analyzed  buckling  of 
radial  spring  supported  rings  and  found  the  critical  ex¬ 
ternal  buckling  pressure  to  be 

P,r- 2(KR)EbI/R3)  [Eq  1 5] 

As  a  first  approximation  in  Ihcsc  analyses  of  cylinders 
and  rings  on  clastic  foundations,  only  bending  defor¬ 
mations  arc  considered,  and  deformations  due  to  cir- 
cumferenlia!  compression  forces  arc  ignored. 

Cylinders  supported  by  a  continuous  clastic  medi¬ 
um:  Hoeg34  gives  an  elastic  solution  for  stresses  and 
displacements  of  a  cylinder  in  an  clastic  medium  con¬ 
sidering  N  -  «j|,/ovas  a  variable  for  the  cases  of  perfect 
bond  and  of  free  slip  between  the  cylinder  and  the 
medium.  He  finds  that  the  behavior  is  a  function  of  the 
relative  compressibility  and  flexibility  of  the  cylinder 
and  the  medium  as  given  by  a  flexibility  ratio 


1  /(»( l-2u/l-u)M 

e,,i/r3ov-) 


[Eq  16] 


and  a  compressibility  ratio 

_ _ M/d-u) 

C  ~  Et/R(lV) 


[Eq  17) 


where  M  is  the  one  dimensional  constrained  modulus 
of  lire  medium  and  v  is  the  Poisson’s  ratio  of  the  medi¬ 
um.  Note  that  the  significant  cylinder  parameters  arc 
the  same  as  those  for  the  unsupported  cylinder. 


' 1  M.  Cheney,  “Bending  anil  Buckling  of  Thin-Walled  Open- 
Section  Rings,”  Proc.  ASCF,  Journ.  I’ng.  Meeh.  Div.  LM5, 
Paper  3665  (1963). 

’ 3  U.  Luscher,  Study  of  the  Collapse  of  Small  Soil-Surrounded 
Tubes,  Technical  Report  Ah'S WC-TDR-6 3-6  (Air  force 
Special  Weapons  Center.  1963). 

54  K.  Hoeg,  Pressure  Distribution  on  Underground  Structural 
Cylinders.  Technical  Report  No.  AEW1.-TR-65-98  (Air 
f  orce  Weapons  Uiboratorv,  1966). 
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Cylinders  in  soil  materials:  On  the  basis  of  empiri¬ 
cal  data  from  tests  of  buried  culverts.  Spangler5  5  devel¬ 
oped  the  well-known  “Iowa  formula”  in  which  the 
detormalions  ot  a  thin-walled  buried  cylinder  due  to 
surcharge  loading  may  be  expressed  in  terms  of  p/Evl / 
R3  and  eR/Eel/R3  where  p  is  a  measure  of  the  sur¬ 
charge  pressure  acting  on  the  cylinder,  and  c  is  a  pas¬ 
sive  deformation  modulus  of  the  soil.  Watkins  and 
Spangler56  present  arguments  indicating  that  the 
quantity  “cR”  and  not  just  “e”  is  a  constant  for  a 
given  soil. 

Luscher57  analyzes  deformation  data  from  tests 
on  buried  cylinders  in  terms  of  (he  parameters  p/Ey  1/ 
R3  and  E'E8I/R3,  where  E'  is  a  deformation  modulus 
of  the  soil  which  he  found  to  correlate  with  M,  the 
one-dimensional  constrained  deformation  modulus. 

It  is  interesting  to  note  that  Spangler's  Iowa  for¬ 
mula  can  be  written  in  a  form  which  is  very  similar  to 
the  first  term  in  the  Hetcnyi  solution  for  the  behavior 
of  rings  on  clastic  foundations. 

Hocg  ran  tests  on  cylinders  buried  in  Ottawa  sand 
which  agreed  well  with  predictions  from  his  elastic 
solution  for  the  behavior  of  a  cylinder  in  an  clastic 
medium. 

Luscher58  summarizes  data  on  the  buckling  of 
cylinders  in  soil  accumulated  by  a  number  of  investi¬ 
gators  and  found  that  the  bending  stiffness  E8I/R3  of 
the  cylinder  correlated  with  failure  by  buckling. 

A  review  of  the  preceding  discussions  shows  that 
the  most  significant  cylinder  parameter  governing  the 
behavior  of  a  cylinder  subjected  to  external  pressures  is 
the  bending  stiffness  E8l/R3(l-u82 ).  in  the  special  case 
in  which  the  external  loading  closely  approaches  a  uni¬ 
form  radial  pressure,  the  hoop  compression  defor¬ 
mation  is  more  important  than  the  bending  defor¬ 
mations,  and  the  important  parameter  becomes 
E8t/R(l-u82).  The  quantity  (l-u82)  is  quite  close  to 
unity  and  may  be  disregarded  with  little  error. 


1  *  MLG.  Spangler,  The  Structural  Design  of  Flexible  Pipe  Cul¬ 
verts,  Bull.  1 53,  Iowa  F.ng  F.xp.  Sla.  (1941). 

**  R.K.  Watkins  and  M.CI.  Spangler,  “Some  Characteristics  of 
the  Modulus  and  Passive  Resistance  of  Soil,  A  Study  in 
Similitude,"  HRB  Proe.  Vol  37  (1958)  pp  576-583. 

* 7  U.  Luscher.  Behavior  of  Flexible  Underground  Cylinders. 
Technical  Report  No.  AI'WL-TR-65-99  (Air  Force  Wcpons 
Laboratory,  1965). 

*  *  U.  Luscher.  Behavior  of  Flexible  Underground  Cylinders. 


Hence,  for  the  purposes  of  dimensional  analysis 
and  modeling,  the  stiffnesses  E0 1/R3  and  E8t/R  arc 
considered  the  significant  parameters  describing  the 
tunnel  lining,  rather  than  the  individual  quantities  E,  I, 
R,  and  t.  This  has  the  great  advantage  of  allowing  the 
model  lining  to  be  made  of  any  material,  regardless  of 
the  prototype  material,  since  it  is  necessary  to  consider 
only  the  structural  stiffness  properties  of  the  liner  and 
not  the  actual  liner  material  properties.  This  assumes 
that  the  prototype  lining  is  subjected  only  to  stresses  in 
the  pseudo-clastic,  working-stress  range,  so  that  the 
strength  parameters  of  the  liner  material  need  not  be 
modeled.  This  is  a  reasonable  assumption  for  most  tun¬ 
nel  linings.  It  must  be  noted,  however,  that  the  ratio 
between  the  two  stiffnesses  is  distorted  somewhat  if 
the  model  lining  is  constructed  of  a  material  different 
from  that  of  the  prototype  lining.  This  means  that  the 
response  of  a  model  lining  will  not  reproduce  with 
complete  accuracy  the  response  of  a  prototype  lining 
constructed  of  a  different  material. 

Variables  describing  the  response  of  the  system. 

The  radial  movements  of  the  tunnel  wall  arc  given  by 
u.  strains  at  points  in  the  rock  mass  behind  the  wall  are 
given  by  c.  and  stresses  within  the  rock  mass  arc  given 
by  a.  The  subscripts  0  and  r  will  be  used  with  these 
variables  to  indicate  the  circumferential  and  radial  di¬ 
rections.  respectively. 

The  response  of  the  tunnel  liner  is  given  by  the 
displacements  u,  and  by  thrusts  T  and  moments  M  in 
the  liner.  The  thrusts  and  moments  arc  given  per  unit 
length  of  tunnel:  pounds  per  inch,  and  inch-pounds  per 
inch.  The  actual  stresses  in  the  liner  arc  not  considered 
because  they  arc  of  less  fundamental  significance  than 
the  thrusts  and  moments.  The  thrusts  and  moments  are 
dependent  only  upon  the  more  general,  or  “first- 
order."  variable  such  as  flexural  and  compression  stiff¬ 
ness.  The  actual  stresses,  however,  arc  also  dependent 
upon  less  significant,  or  “second-order,”  details  of  the 
lining  design. 

Dimensionless  Pi  Terms.  The  behavior  of  the  proto 
type  tunnel  is  assumed  to  be  determined  and  described 
with  sufficient  accuracy  by  the  variables  giver  in  Table 
A2.  These  variables  are  the  Xj  terms  of  Equation  1, 
which  then  becomes 

'(»V,  °h,  c  or  q„,  <p ,  ot,E,  u.  d  or  a  or  D  or  R. 

0j.  (7,  s.  Ejl/R* ,  E8t/R,u.  c,  0,  T,  M)  =  0 

|Eq  18) 
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Table  A2 

Dimensionless  Pi  Terms 

Loading: 

°h /°v  °vAlu 

Intact  rock: 

0.  I'-Ali,.  i> 

Rock  mass: 

0j.  0,  D/s 

Support  properties:  KRJ/1-Vt  or  «VR3/Iivl.  KR/I-m  or  ovR/L,,| 

Response : 

u/D.  t .  o/ov  or  <J/n„,  T/<»VR,  M/ovR5 

The  next  step  in  the  determination  of  the  simili¬ 
tude  requirements  governing  a  model  study  of  this 
phenomenon  is  to  perform  a  dimensional  analysis;  that 
is.  to  determine  u  set  of  pi  terms  as  in  Equation  2.  By 
inspection  of  Table  A1  and  Equation  3  it  is  seen  that 
there  arc  m  =  n  -  r  =  18  -  2  -  16  dimensionless  pi 
terms  which  describe  the  phenomenon.  Rigorous  meth¬ 
ods  arc  available  for  determining  a  complete  set  of 
independent  pi  terms  (see  for  example,  Langhaar  or 
Murphy).  However,  a  complete  set  for  the  variables  in 
Table  A I  can  be  determined  by  inspection,  and  arc 
given  in  Table  A2. 

Note  that  in  several  places  in  Table  A2  two  ex¬ 
pressions  arc  given  as  being  equally  valid.  The  princi¬ 
ples  of  dimensional  analysis  state  that  in  any  given  set. 
all  pi  terms  are  independent;  i.c.,  no  one  pi  term  can  be 
obtained  by  multiplying  together  any  combination  of 
the  other  pi  terms  of  that  set.  But  there  are  an  infinite 
number  of  sets  of  independent  pi  terms  which  can  be 
given  for  the  n  variables,  any  one  set  of  which  can  be 
derived  from  any  second  set  by  multiplying  terms  of 
the  second  set  together.  For  examples,  either  o/ov  or 
a/qu  may  be  considered  as  one  of  the  pi  terms,  but 
both  cannot  be  considered  as  part  of  one  set  because 
they  are  not  independent.  One  can  be  derived  from  the 
other  by  multiplication  with  other  of  the  pi  terms  as 
follows: 

o/ov  =  o/qu  x  qu/av  |Eq  19] 

It  was  mentioned  earlier  that  much  insight  into  the 
phenomenon  being  considered  could  often  be  obtained 
from  the  relationships  of  the  variables  in  the  pi  terms. 
The  terms  ov/qu  and  0|,/ov  show  that  the  frcc-field 
stress  level  relative  to  the  strength  of  the  material  is 
important,  and  that  the  ratio  of  horizontal  to  vertical 
stress  is  of  significance.  The  term  D/s  shows  that  the 
ratio  of  tunnel  diameter  to  joint  spacing  is  significant 
in  the  behavior  of  the  tunnel.  The  ratios  ER3/EVI  and 
ER/E«t  (E/Evl/R3  and  E/Evt/R)  show  that  die  relative 


stiffnesses  of  the  rock  mass  and  the  liner  arc  of  impor¬ 
tance  in  influencing  the  behavior.  Alternately,  from 
ovR3/Eel  and  ovR/Evt  it  is  seen  that  the  behavior  of 
the  liner  depends  upon  its  stiffness  relative  to  the  ap¬ 
plied  stress.  The  development  of  thrusts  and  moments 
is  best  expressed  in  terms  of  the  parameters  T/ovR  and 
M/ovR2.  The  importance  of  many  of  these  terms,  of 
course  has  been  recognized  for  some  time  from  other 
work,  but  this  dimensional  analysis  is  an  independent 
substantiation  of  this  importance. 

Development  of  Modeling  Laws.  Using  die  pi  terms 
just  derived.  Equation  2  becomes 

F(r;h/ov,  r;v  /c)u.  ot/qu.  0,  E/q„,  u,  ER3/EtI 
or  ovR3/E(jl.0j,  0,  D/s.  ER/Eot 

1  Eq  201 

or  rjvR/Evt,u/d,  c,  o/av 
or  o/qu.  T/ovR,  M/ovR2)  -■  0 

This  is  the  dimensionless  functional  relationship  de¬ 
scribing  the  phenomenon.  If  the  model  study  is  to 
accurately  reproduce  the  prototype  field  behavior,  it  is 
necessary  that  the  pi  terms  as  given  it'.  Table  A2  and 
Equation  20  be  identical  for  the  model  and  the  proto¬ 
type.  as  indicated  by  Equation  4.  This  is,  for  example: 

(O|,/ov)  mode!  =  (cj|,/<jv)  prototype  ^  ^ 

(0)  mode!  =  (0)  prototype 
and  so  on  for  the  rest  of  the  pi  terms. 

If  the  symbol  KXj  is  used  to  represent  the  ratio 
between  the  value  of  one  of  the  xj  terms  in  the  model 
and  in  the  prototype,  for  example: 

Kov  =  (ov)  tnodel/(ov)  prototype  [Eq  22] 

then  the  requirements  of  similitude  as  given  in  Equa¬ 
tions  4  and  21  dictate  certain  relationships  which  must 
exist  between  the  KXj  ratios.  The  KXj  ratios  are  called 
the  scale  factors,  and  the  relationships  between  them 
arc  called  the  model  laws.  The  model  laws  for  this 
phenomenon,  as  derived  from  the  pi  terms  of  Table  A2 
by  simple  algebraic  manipulation,  are  given  in  Table 
A3. 

The  .node!  laws  of  Table  A3  show  that  all  of  the  Xj 
quantities  having  the  dimensions  of  length  L  scale  in 
the  same  ratio  Kg  between  the  model  and  the  proto¬ 
type:  all  Xj  quantities  having  the  dimensions  of  stress 
FL’2  scale  in  the  same  ratio  K0;  and  all  dimensionless 
quantities  such  as  strains  and  angles  have  the  same 
magnitude  in  the  model  as  in  the  prototype.  Further¬ 
more.  the  model  laws  show  that  the  scale  factors  for 
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Table  A3 

Model  Laws 

1  invar  Dimensions  t.. 

K|.  *  Kt,  ••  K„  •••  K|>  -  K|<  = 

St  rossvs : 

“  K<7v  -  k0|,  “  Kiju  -  Kat  “  K|.* 

KKwI/K*  =  KkV!/r 

Strains  and  Angles: 

~  k,  -  Kj,  -  K0  r:  l 

K  T  a  K«Ki,,  Km  =  K0KLJ 

Response: 

lengths  and  stresses.  Kl  and  Kf;,  are  independent  and 
may  be  chosen  arbitrarily. 

The  inter-relationships  between  the  scale  factors  Kp. 
K%j,  Ko<  and  K  p  are  not  to  be  interpreted  as  restric¬ 
tions  upon  this  independency  of  K|  and  Kp  but  mere¬ 
ly  as  a  definition  of  the  manner  in  which  thrusts  and 
moments  scale  between  model  and  prototype. 

If  the  model  laws  of  Table  A3  arc  satisfied,  then 
the  requirements  of  Equation  4  arc  satisfied  and  the 
model  is  “similar”  to  the  prototype  and  behaves  exact¬ 
ly  as  the  prototype.  By  measuring  the  \j  quantities  in 
the  model,  the  Xj  quantities  of  the  prototype  can  be 
predicted  through  the  scale  factors.  The  accuracy  with 
which  the  predicted  prototype  behavior  matches  the 
actual  prototype  behavior  depends  upon  two  factors: 

(1)  The  accuracy  of  the  assumption  that  the  Xj 
quantities  of  Table  A1  arc  the  quantities 
which  determine  and  describe  the  phenome¬ 
non,  and 

(2)  The  accuracy  with  which  the  model  laws  arc 
satisfied. 

Model  Rock  Material  Requirements.  Some  very 
stringent  limitations  on  the  behavior  of  the  model  rock 
material  arc  implicit  in  the  modeling  ratios  given  above. 
The  model  laws  K0(  =  Kqu  =  Kay  =  K0h  and  =  I 
require  that  on  any  dimensionless  plot  of  strengths,  the 
data  for  both  the  model  and  prototype  rock  materials 
must  collapse  onto  a  single  curve.  For  example,  if  Fig. 
AS  were  the  Mohr  envelopes  for  the  prototype  and 
model  materials,  then  on  a  dimensionless  plot  of  r/qu 
vs  o/qu  such  as  Fig.  A6  (lie  envelopes  for  the  two 
materials  must  coincide.  On  any  other  dimensionless 
strength  plot,  such  as  (a,  -o3)/qu  vs  o.i/qu-  the 
strength  envelopes  for  the  two  materials  must  also  co¬ 
incide. 

A  similar  relationship  exists  for  the  deformation 
characteristics  of  the  model  and  prototype  materials. 


The  modeling  ratios  K1:  -  Kt|u,  Kt  =  1,  Ky  =  I,  and  the 
strength  requirements  given  above  require  that  the 
materials  have  the  same  Poisson’s  ratio  and  that  on  all 
dimensionless  plots  of  strains  or  deformations  versus 
stress  the  curves  for  (he  two  materials  must  coincide. 
For  example,  il  Fig.  A7  were  the  stress-strain  curves  for 
the  model  and  prototype  materials  at  comparable  con¬ 
fining  pressures;  i.e..  at  a  value  of  Oj/qu  which  is  the 
same  for  both  materials,  then  on  a  dimensionless  plot 
of  (a,  o3 )/ q„  vs  c,  the  curves  for  both  materials  must 
coincide  as  in  Fig.  AS.  This  means,  for  example,  that  in 
triaxial  compression  tests  at  comparable  confining  pres¬ 
sures.  the  materials  must  fail  at  the  same  strains. 
Hence,  on  a  dimensionless  plot  of  c-faj[urcvs  a3/ qu,  the 
data  for  the  two  materials  must  collapse  onto  a  single 
curve. 

In  practice  these  requirements  are  almost  impos¬ 
sible  to  «”tisfy.  In  a  model  study  of  a  linearly  elastic 
phenomenon,  the  strength  modeling  laws  do  not  exist, 
and  the  deformation  modeling  laws  arc  not  so  critical 
since  the  stress-deformaiion  relationships  arc  linear. 
For  example,  K0  =  2  might  bo  dlowcd  without  serious¬ 
ly  affecting  the  accuracy  of  the  model.  For  a  model 
study  in  which  inelastic  deformations  and  failure  con¬ 
ditions  are  important,  however,  the  model  laws  must 
be  satisfied  as  nearly  as  possible.  Patterns  of  stress  and 
strain  distribution  in  the  prototype  may  change  mark¬ 
edly  as  non-linear,  inelastic  deformations  occur  and  as 
failure  conditions  arc  approached  or  reached.  If  the 
model  laws  arc  not  satisfied  and  the  model  materials  do 
not  fulfill  the  requirements  outlined  above,  the  pat¬ 
terns  of  stress  and  strain  distribution  in  the  model  may 
differ  considerably  from  those  of  the  prototype. 

In  modeling  studies  the  most  common  solutions  of 
this  dilemma  are  either  ( I)  to  conduct  tests  in  only  the 
quasi-elastic,  working-stress  range,  where  the  model 
material  requirements  arc  not  so  critical;  or  (2)  to  build 
the  model  from  essentially  the  same  material  as  the 
prototype,  as  is  done  in  microconcretc  model  studies 
of  reinforced  concrete  structures.  These  approaches  arc 
not  feasible  in  geomechanical  model  studies  of  under¬ 
ground  openings,  however. 

In  the  first  case,  underground  structures  arc  highly 
indeterminate  so  that  local  failures  can  develop  with¬ 
out  leading  to  complete  failure  of  the  structure.  In 
fact,  the  economics  of  underground  construction  often 
demand  that  such  local  failures  be  tolerated.  Hence,  a 
study  of  the  low  stress,  quasi-elastic  behavior  is  not 
sufficient. 
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Figure  A5.  Prototype  and  model  Mohr  envelopes. 


Figure  A6.  Dimensionless  prototype  and  model  Mohr  envelopes. 
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Figure  A7.  Prototype  and  model  stress-strain  curves. 


Figure  A8.  Dimensionless  prototype  and  model  stress- 
strain  curve. 


The  second  approach  is  not  feasible  for  three  rea¬ 
sons.  First,  discontinuities  in  the  actual  rock  mass  exist 
on  a  physical  scale  such  that  it  is  impossible  to  obtain 
samples  of  the  rock  mass  small  enough  so  that  accurate 
model  studies  can  he  performed  when  the  effect  of 
mass  discontinuities  is  being  studied.  Secondly,  the 


strength  of  the  rock  materials  is  generally  so  great  that 
the  si/.c  of  models  constructed  of  the  prototype ma¬ 
terial  which  can  economically  be  loaded  to  failure  in 
the  laboratory  are  loo  small  to  be  of  interest.  Finally, 
it  would  be  very  difficult  to  gel  “identical  samples”  for 
testing  from  a  real  rock  mass  because  most  rock  for¬ 
mations  are  not  sufficiently  homogeneous.  Hence  it  is 
necessary  to  use  artificial,  low  strength  materials  for 
the  construction  of  gcnmcchauical  models. 

Data  for  intact  rock  (tor  example,  Deere  and  Mil¬ 
ler.59  Deere/10  Hendron/'  Handin  and  Hager/’2 
Handin  et  al/ '  Corps  of  Engineers/4  and  Robert¬ 
son65  indicate  that  on  the  average,  the  properties  of 
intact  rock  are  such  that  the  tensile  strength  is  about 
five  percent  to  ten  percent  of  the  unconfined  compres¬ 
sive  strength  and  the  modulus  of  elasticity  is  about  250 
to  500  limes  the  unconfined  compressive  strength  (E 
being  the  tangent  modulus  at  fifty  percent  of  <|u), 
while  the  angle  of  internal  friction  commonly  varies 
between  25°  and  60°,  and  v  is  between  0.1  and  0.3. 
That  is,  for  actual  rock: 

5'a<oi/i|u<  HY/r :  200  <  E/q„  <  500; 

25°  <  <>  <  60" ;  0. 1  <  v  <  0.3 

The  material  chosen  for  the  construction  of  the  model 
should  also  have  properties  within  these  ranges  if  simili¬ 
tude  is  to  be  achieved. 

Test  data  indicate  that  rock  specimens  typically 


*D.V.  Deere  and  H  P.  Miller,  Engineering  Classification  and 
Index  Properties  for  Intact  Rock,  Technical  Report  No. 
Af  WL-TR-65-1  !(>  (Air  force  Weapons  Laboratory.  1966). 

*  "  D.V  Deere,  “Geologic  Considerations.” 

-  1  A..I.  I  lend  roii,  “Mechanical  Properties  of  Rock.” 

‘ 5  J.  Handin  and  R.V.  I  lager,  “Experimental  Deformation  of 
Sedimentary  Rocks  Under  Confining  Pressures:  Tests  at 
Room  Temperature  on  Dry  Samples."  Ik  ill  A.A.P.G.,  Vol 
41.  No.  I  (January  1957)  pp  1-50. 

*JJ.  Handin.  R.V.  Hager,  M.  friedinan  and  J.N.  feather,  “Ex¬ 
perimental  Deformation  of  Sedimentary  Rocks  Under  Con¬ 
fining  Pressure:  Pore  Pressure  Results,”  Bull.  A.A.P.G.,  Vol 
47.  No.  5  (May  1963)  pp  717-755. 

*  *  Tests  for  Strength  Characteristics  of  Rock,  Pile  Driver  Proj¬ 

ect  MRD  l.al>  No.  64/90  (U.S.  Army.  Missouri  River  Divi¬ 
sion  Laboratory  |MRD|.  Sept  1964);  Tests  for  Strength 
Characteristics  of  a  Schistose  Gneiss,  MRD  Lab  No.  64/493 
(May  I9(>5);  Strength  Parameters  of  Selected  Intermediate 
(duality  Rocks,  MUD  Lab  No.  64/493  (July  1966). 

'  '  f.C.  Robertson,  "I  xpeiiment.il  Study  of  the  Strength  of 
Rocks,  Hull.  C.S.A..  Vol  66  (October  1955)  pp  1275-1314. 
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fail  in  unconfincd  compression  at  axial  strains  of  0.2  to 
1.0  percent.  In  triaxial  compression  at  confining  pres¬ 
sures  equal  to  their  unconfined  compression  strength, 
rock  specimens  reach  a  peak  stress  difference  (oi-oj) 
at  axial  strains  which  may  range  widely,  from  around 
one  percent  for  dense  igneous  rocks  up  to  ten  percent 
to  twenty  percent  or  more  for  ductile  shales  or  evapo- 
rites.  Rocks  commonly  exhibit  dilation  during  shear, 
possible  exceptions  being  very  porous  sedimentary  or 
volcanic  rocks  whose  porous  structure  collapses  during 


shear,  or  some  evaporites  which  may  fail  by  inter- 
crystalline  gliding  with  no  volume  change.  A  “general 
rock  modeling  material”  should  exhibit  dilation  during 
shear  to  satisfy  similitude. 

Typical  stress-strain  curves  for  a  wide  range  of 
rock  types  are  given  in  the  references  listed  above.  The 
modeling  material  chosen  should  possess  stress-strain 
curves  which  arc  of  the  same  shape  as  those  of  a  typical 
rock,  as  in  Fig.  A4. 
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APPENDIX  B:  DESCRIPTION  OF  MODEL  LOADING  APPARATUS* 


Lateral  Loading  Elements.  The  design  criteria  for 
tlie  system  ol  applying  the  lateral  pressures.  r»v  and 
to  the  model  were  that  it  be  able  to  applv  a  uniform 
pressure  up  to  1.000  psi  on  the  24"  x  X"  laces  ol  the 
model  lor  a  total  reaction  of  l')2.00l)  pounds,  and  that 
it  be  able  to  deform  up  to  one-ipiarler  inch  as  the 
model  was  loaded.  Although  large,  irregular  defor¬ 
mations  of  the  faces  were  not  anticipated,  a  nonuni¬ 
form  concave  outward  deformation  was  anticipated  for 
two  reasons:  (I)  the  presence  of  the  tunnel,  and  (2) 
restraining  friction  along  the  adjacent  loading  faces. 
Hence,  the  lateral  brails  had  to  be  applied  hv  a  flexible 
leading  system  which  could  adjust  to  the  irregular 
model  deformations,  rather  than  by  a  rigid  loading 
head. 

It  was  decided  that  the  load  should  be  actively 
applied  to  all  four  lateral  sides,  rather  than  actively 
loading  the  model  on  two  adjacent  sides  and  pushing  it 
against  a  passive  reaction  on  the  opposing  sides.  Al¬ 
though  this  increased  the  complexity  of  the  loading 
system,  it  was  done  to  maintain  loading  symmetry  in 
the  model  so  that  friction  between  the  model  and  the 
longitudinal  and  lateral  loading  heads  would  be  sym¬ 
metrical  about  the  tunnel. 

The  loading  system  which  was  chosen  is  a  com¬ 
pletely  mechanical  true  whose  basic  operation  is  quite 
similar  to  a  system  used  by  Hock66  on.  small  models  6” 
x  6"  x  I".  It  corrsists  of  a  pyramid  of  increasingly 
larger  triangular  elements  produced  by  welding  angle 
irons  and  Hat  plates  together,  as  shown  in  Fig.  111.  Two 
such  sets  of  elements  arc  used  to  apply  the  load  to  each 
lateral  face  of  the  model,  as  shown  in  Fig.  B2.  The  load 
is  applied  by  a  hydraulic  jack  against  the  transition 
head  clement  and  distributed  down  through  the  pyra¬ 
mid  to  the  model,  which  is  in  contact  with  the  smallest 
triangular  elements,  Element  No.  4.  At  the  contact  be¬ 
tween  the  eleme  grooves  I / 1 f* -in.  deep  with  a  I -in. 
radius  arc  cut  into  the  plates  and  the  angle  corners  arc 
rounded  on  a  l/l-in.  radius.  Tire  purpose  of  the 
grooves  and  rounded  corners  is  threefold:  (I)  to  allow 
the  elements  to  rotate  with  respect  to  each  other  and 
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to  adjust  to  any  irregular  deformations  of  the  model 
while  carrying  equal  loads.  (2)  to  aid  in  aligning  the 
elements,  and  (3)  io  provide  a  large  enough  bearing 
area  at  the  contacts  to  prevent  local  yielding  of  the 
steel  and  flattening  of  the  contacts.  Even  so,  high  con¬ 
tact  stresses  exist,  and  a  high  yield  strength  steel  (T-l 
steel)  was  used  for  the  transition  heads  and  the  angles 
of  Element  No.  1.  This  system  of  applying  the  lateral 
loads  was  chosen  for  two  reasons:  (1)  the  successful 
operation  of  the  similar  loading  system  used  by 
Hock/'7  and  (2)  the  simple,  completely  mechanical 
nature  of  Hie  system,  which  suggested  that  it  would  be 
very  dependable  and  tugged. 

Load  Distribution  Characteristics  of  the  Triangular  Ele¬ 
ment  Lateral  Loading  System.  The  load  distribution 
characteristics  of  the  triangular  loading  clement  as¬ 
sembly  were  studied  in  some  detail  to  determine  how 
well  they  satisfy  the  lateral  boundary  condition;  that 
of  a  uniform  lateral  stress  distribution  some  distance 
from  the  tunnel.  It  is  recognized  that  the  stress  distri¬ 
bution  on  the  face  of  the  model  will  be  quite  irregular 
due  to  the  finite  width  and  the  finite  stiffness  of  the 
loading  elements.  The  philosophy  guiding  the  design  of 
tire  dements  was  that  they  should  be  relatively  narrow, 
stiff,  and  closely  spaced,  and  that  cadi  should  apply 
the  same  total  force  to  the  model  even  if  it  deforms 
unevenly.  Then,  within  a  short  depth  into  the  model, 
approximately  equal  to  the  width  >f  the  elements,  the 
actual  stress  distribution  would  deviate  only  slightly 
from  the  average  stress  applied  to  the  boundary.  To 
check  the  behavior  of  the  loading  apparatus,  three 
things  were  done.  First,  one  set  of  loading  elements 
was  instrumented  to  determine  the  loads  being  carried 
by  the  smallest  triangles,  under  uneven  deformations  of 
the  assembly.  Second,  concrete  blocks  the  size  of  the 
actual  model  blocks,  24”  x  24”  x  8”.  both  with  and 
without  tunnels,  were  tested  m  biaxial  compression. 
These  were  sprayed  with  a  brittle  lacquer  coating  on  an 
unloaded  24”  x  24”  face  to  study  the  strain  distri¬ 
bution  in  the  block.  Third,  before  testing  the  first 
model  tunnel,  ,t  solid  block  of  the  model  material  with¬ 
out  a  tunnel  was  tested  with  internal  instrumentation 
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Pyramid  of  Triangular  Clements, 
7  Per  s;dr 


Figure  B2.  Sketch  of  lateral  loading  element  assembly  (not  to  scale). 


to  measure  the  strain  field  which  was  produced  within 
the  block.  These  studies  led  to  the  conclusion  that  the 
lateral  loading  system  was  performing  satisfactorily  and 
fulfilling  the  design  requirements  and  boundary  con¬ 
dition. 

Load  Distribution  As  Measured  by  the  Instrumented 
Elements.  One  set  of  triangular  elements  as  shown  in 
Figure  B4  was  used  for  this  study.  An  electrical  resist¬ 
ance  strain  gage  (SR4  Gage,  Type  A-7)was  placed  on  the 
center  of  each  leg  of  each  of  the  four  small  triangular 
elements  (Element  No.  3)  as  shown  in  Fig.  B3.  The  two 
gages  from  each  clement  were  wired  into  opposite  arms 
of  a  four-arm  Wheatstone  bridge  and  monitored  by  a 
strain  indicator.  Each  of  the  small  triangular  elements 
was  then  loaded  individually  with  the  transition  head 
to  obtain  a  calibration  curve  showing  SR4  gage  reading 
versus  total  load  carried  by  the  individual  clement.  The 


entire  pyramid  of  triangular  loading  elements  and  tran¬ 
sition  head  was  then  loaded  and  at  successive  stages  of 
loading  the  strain  indicated  by  the  SR4  gages  on  each 
of  the  small  elements  (No.  3)  was  recorded  and  com¬ 
pared  with  the  calibration  curves  to  determine  how 
much  of  the  total  applied  load  was  being  carried  by 
each  individual  clement. 

This  method  of  monitoring  the  loads  carried  by 
individual  elements  is  not  satisfactory  for  general 
usage,  although  it  was  sufficient  for  the  purpose  of  the 
immediate  investigation.  The  problem  with  this  meth¬ 
od  is  that  it  is  very  sensitive  to  bending  moments  in  the 
triangle  legs,  which  are  in  turn  caused  by  bending  of 
the  bottom  flat  plate  of  the  clement  at  the  contact 
with  the  model.  This  system  of  instrumentation  is  thus 
quite  sensitive  to  pressure  distribution  on  the  base  of 
the  elements.  Ideally,  the  instrumentation  should  be 
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sensitive  only  in  tin*  <t\ei  ..  . . n  t i is- 

legs  of  the  angles.  flu-.  imgln  !  ..  ,.r  s.ifiole. 

hy  placing  another  A- 7  ga  nv  ...  i. -v  Iv  h  .mi!|0 
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For  the  immediate  investigation,  this  difficulty 
was  solved  as  shown  in  Fig.  B4.  A  steel  bar  j/2”  x 
i/s"  was  placed  under  each  edge  of  each  clem... it  and 
these  m  turn  were  placed  upon  a  I /S  in.  strip  of  rub¬ 
ber.  Rubber  snips  of  different  stiffness  were,  used  to 
induce  differentia!  movement  of  the  elements  during 
loading.  The  use  of  narrow  strips  under  the  edges  of 
the  elements  eliminated  the  problem  of  pressure  distri¬ 
bution  across  the  base.  In  addition  the  elements  were 
placed  upon  a  12"  x  8”  x  1-1/2”  steel  plate  which 
could  be  moved  freely  so  that  the  individual  elements 
could  be  centered  in  the  test  machine  for  calibration, 
and  then  the  whole  assembly  could  be  centered  for 
loading.  In  this  way  the  whole  loading  assembly  could 
be  tested  with  the  small  elements  in  exactly  the  same 
position  and  with  exactly  the  same  base  conditions 
with  which  they  were  calibrated. 

The  results  of  the  tests  on  the  whole  assembly  are 
shown  in  Fig.  B5.  The  maximum  variation  between 
elements  is  approximately  plus  or  minus  2.2  percent  of 
the  total  assembly  load  and  was  generally  around  plus 
or  minus  one  percent.  The  variation  between  the  load 
carried  by  each  element  is  almost  equal  to  the  variation 
in  the  calibration  of  the  individual  elements,  so  that 
within  the  accuracy  of  this  method  of  mca -urement, 
the  elements  are  carrying  equal  loads.  This  was  meas¬ 
ured  during  tests  in  which  significant  movement  of  the 
whole  assembly  and  significant  relative  movements  of 
individual  elements  were  recorded  as  shown  in  Fig.  B5. 

In  summary,  this  study  showed  that  this  lateral 
loading  assembly  can  accomodate  substantial  differ¬ 
ential  movements  of  its  base  while  distributing  the 
total  load  equally  across  the  base. 

Longitudinal  Restrain  and  Reaction  Frame.  The  meth¬ 
od  of  applying  the  longitudinal  restraint  and  the  meth¬ 
od  of  supplying  the  lateral  jacking  reactions  arc  inti¬ 
mately  related  in  the  design  of  the  loading  machine  and 
must  be  considered  together.  As  discussed  previously,  a 
condition  of  plane  strain  should  be  approximated  in 
the  model.  Three  methods  of  achieving  this  result  were 
listed: 

(1)  A  uniform  pressure  against  the  longitudinal 
faces,  controlled  to  null  any  longitudinal 
expansion  which  tends  to  develop. 

(  2)  Rigid  heads  against  the  longitudinal  faces,  tied 
rigidly  together  across  the  model. 

(2)  Rigid  heads  against  the  longitudinal  faces, 
with  a  controlled  load  applied  to  them  to  null 


any  longitudinal  expansion  which  tends  to 
develop. 

Uniform  Restraining  Pressure  System  It  is  believed 
that  the  uniform  pressure  method  would  not  satisfy 
the  boundary  condition  satisfactorily.  Dilation  of  the 
material  is  expected  in  failure  /ones  which  develop 
around  the  tunnel,  with  some  expansion  occurring  in 
the  longitudinal  direction.  More  pressing  would  be  re¬ 
quired  to  null  longitudinal  strains  in  these  regions  than 
in  regions  away  from  the  tunnel  which  are  still  behav¬ 
ing  elastically.  A  uniform  nulling  pressure  would  then 
be  too  low  to  prevent  expansion  in  the  plastic  /ones 
and/or  too  high  to  exactly  null  strains  in  the  elastic 
zones.  Hence,  the  longitudinal  faces  would  not  remain 
plane  during  the  loading,  but  would  warp.  Thus,  (he 
system  is  unsatisfactory. 

Rigid  Heads  Tied  Across  Model.  The  second  method 
suggested  has  the  obvious  advantage  ol  greatly  simpli¬ 
fying  the  test  procedure.  Once  the  model  is  in  testing 
position,  the  longitudin  d  heads  are  brought  into  con¬ 
tact  with  it  and  tied  together  across  the  model,  prob¬ 
ably  at  the  corners  of  the  model.  Then  during  the 
lateral  loading  the  model  is  restrained  from  longitudi¬ 
nal  expansion  by  the  rigid  heads  and  ties.  This  would 
eliminate  the  need  for  monitoring  and  regulating  the 
longitudinal  deformations  and  loads.  A  serious  problem 
which  must  be  given  careful  consideration  with  an  ap¬ 
paratus  of  this  sort  is  the  difficulty  of  seating  with  the 
restraining  heads  against  the  model.  Because  of  the 
very  small  longitudinal  expansion  which  would  occur 
even  in  the  plane  stress  condition,  the  restraining  heads 


must  he  seated  very  carefully  ami  lightly.  Otherwise, 
the  expansion  of  the  model  which  would  occur  before 
intimate  contact  with  the  restraining  loads  was  de¬ 
veloped  would  he  of  the  same  order  of  magnitude  as 
the  total  expansion  during  the  test.  In  such  a  ease  the 
actual  test  condition  would  approximate  plane  stress 
more  closely  than  plane  strain.  In  spite  of  this  problem, 
the  anticipated  simplicity  of  the  testing  procedure  with 
an  apparatus  of  this  sort  is  very  attractive  ami  the  de¬ 
sign  of  such  a  system  was  carefully  considered. 

To  be  acceptable,  the  calculated  deformation  of 
such  an  apparatus  should  not  exceed  about  10  percent 
of  the  plane  stress  expansion  of  a  model  block.  Such  an 
apparatus  proved  too  massive  to  be  practical.  The  prob¬ 
lem  of  adequately  seating  a  model  in  an  apparatus  of 
this  design  remains  unsettled.  Furthermore,  such  an 
apparatus,  representing  a  sizable  investment,  would  not 
have  die  adaptability  and  general  capabilities  of  the 
design  finally  chosen.  For  these  reasons,  it  was  decided 
that  rigid  heads  tied  together  across  the  model  were 
not  a  satisfactory  design. 

Controlled  Rigid  Longitudinal  Heads.  The  next  system 
considered  is  shown  schematically  in  Fig.  Bb.  The  bot¬ 
tom  frame  and  lateral  reactions  would  be  combined 
into  a  single  unit.  The  top  hear!  would  consist  of  a  grill 
wotk  as  shown  in  Fig.  B7  with  a  I -in.  bottom  cover 
plate  in  contact  with  the  model  and  a  1/2-in.  top  cover 
plate.  Hydraulic  jacks  would  force  the  head  against  the 
model  to  null  longitudinal  expansion  of  the  model. 

Calculation  of  the  deformation  of  this  “rigid” 
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Figure  H6  Sketch  of  possible  cantilevered  l.iteial  leactiuti  sssteni  (not  to  scale). 
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Section  AA 

Figure  B7.  Plan  and  section  o|  top  head.  Window  is  I”  x  S”  x  X"  plate  ol  LOI-  I ul lex  glass.  In  plan  view,  crosses 
mark  centers  of  hydraulic  jacks,  circles  represent  i/j"  holes  m  top  plate  ter  ie:eive  1.4  "  rod  threaded 
into  bottom  covei  plate 


r 


!k\kI,  assuming  a  unironnly  distributed  load  from  the 
model  and  point  loads  from  the  jacks,  indicates  that 
waviness  of  the  surface  in  contact  with  the  model 
would  he  less  than  plus  or  minus  It)  percent  of  the 
longitudinal  model  expansion  which  would  occur  in 
the  plane  stress  case,  considering  both  shearing  and 
bending  deforniadnus.  This  was  considered  satisfactory 
and  was  actuary  an  overestimate  since  the  jack  re¬ 
actions  would  be  spread  over  the  area  of  the  jack  bases 
and  would  not  be  point  loads  . 

The  difficulty  with  this  design  'as  that  the  sup¬ 
port  for  the  lateral  jacks  was  not  symmetrical,  but  was 
cantilevered  off  the  bottom  frame.  This  produced 
bending  stresses  and  a  convex  upward  curvature  in  the 
bottom  frame  under  the  base  of  the  model,  which  was 
intolerably  large.  A  number  of  schemes  were  con¬ 
sidered  for  reducing  this  curvature  but  calculations  for 
all  schemes  showed  that  the  maximum  deformation  of 


reactions  were  symmetrically  supported  and  defor¬ 
mation  ot  the  reaction  Irame  was  not  a  major  problem. 
This  design  is  shown  in  Figures  B8  and  B14.  Overall 
views  of  the  apparatus  are  shown  in  Figure  B8  and  B9. 

I  he  lateral  loads  w'ere  applied  by  the  pyramids  of  tri¬ 
angular  clerne  its  discussed  previously,  which  were  ar¬ 
ranged  as  shown  in  Fig.  BIO.  The  reactions  for  the 
lateral  jacks  were  supplied  as  shown  in  Fig.  Bll. 
The  horizontal  1  3/4-in.  lods  were  removed  while  posi¬ 
tioning  the  model  in  the  loading  frame  and  setting  up 
the  test.  These  rods  had  the  same  cross-sectional  area  as 
the  horizontal  3”  x  3"  x  7/1  o  ”  angle  iron  ties,  and  the 
axis  of  the  lateral  jacks  was  spaced  vertically  midway 
between  the  centroid  of  the  rods  and  the  centroid  of 
the  angles,  in  both  the  ov  and  0|,  directions.  Hence, 
bending  moments  in  the  lateral  reaction  frame  were 
kept  to  a  minimum.  Moreover,  the  lateral  reaction 
system  was  independent  of  the  longitudinal  system, 
and  deformations  and  moments  in  the  lateral  system 
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Figure  B12.  Lower  longitudinal  reaction  head  and  con¬ 
crete  pedestal. 

top  and  bottom  longitudinal  reaction  heads  as  shown 
in  Fig.  B14,  which  were  tied  together  by  the  tour  verti¬ 
cal  1  3/4-inch  diameter  rods.  While  considerable  defor¬ 
mation  of  these  top  and  bottom  reaction  heads  oc¬ 
curred,  it  was  of  no  significance  because  the  model  did 
not  rest  directly  against  them.  The  model  was  shielded 
from  irregular  deformation  of  the  bottom  frame  by  the 
2-ft  concrete  cube,  while  being  loaded  from  the  top  by 
the  hydraulic  jacks  acting  against  the  top  loading  head. 
The  system  designed  for  monitoring  the  longitudinal 
strains  is  discussed  later  in  the  section  on  instru¬ 
mentation  of  the  model. 

For  those  members  in  which  deformation  was  not 
a  controlling  factor,  the  reaction  frame  was  designed 
for  an  extreme  fiber  stress  of  20,000  psi  in  both  ten¬ 
sion  and  compression  in  the  rolled  steel  sections,  which 
were  of  A36  steel.  The  horizontal  and  vertical  tic  rods 


were  subjected  to  higher  stresses  in  the  threaded  sec¬ 
tions,  and  were  of  higher  yield  steel. 

In  an  attempt  to  ensure  seating  of  the  loading 
heads  (and  longitudinal  deformation  measuring  points) 
against  the  model,  the  following  procedure  was  used: 
The  top  3/8-in.  aluminum  plate  beneath  the  model 
(Fig.  B13)  was  the  bottom  of  the  mold  in  which  the 
model  was  compacted,  hence  it  fitted  tightly  against 
the  model.  The  bottom  3/8-in.  steel  plate  was  cast  di¬ 
rectly  against  the  concrete  cube.  As  the  model  was  put 
in  testing  position  a  thin  film  of  plaster  (approximately 
l/16-in.  thick)  was  placed  between  the  3/8-in.  plates  to 
ensure  complete  contact  between  them.  Another  thin 
film  of  piaster  was  placed  between  the  top  loading 
head  and  the  model  to  ensure  a  tight,  continuous  fit. 
Then,  when  testing  was  begun,  a  30  psi  longitudinal 
seating  load  was  applied  before  the  lateral  loads  were 
applied. 

The  hydraulic  jacks  used  were  twelve  Simplex 
RC-6010  double-ai  ting  60  ton  hydraulic  rams.  They 
were  actuated  by  the  pressure  console  seen  in  Fig.  B8. 
The  console  was  driven  by  air  pressure  and  featured 
two  independent  hydraulic  systems  each  capable  of 
producing  pressures  of  10,000  psi  from  an  air  pressure 
of  100  psi.  The  four  jacks  applying  the  ov  loads  to  the 
model  were  driven  by  one  of  the  consoles’  hydraulic 
systems,  while  the  four  applying  the  oh  loads  were 
driven  by  the  other  system.  A  hydraulic  pressure  of 
7.150  psi  was  required  to  develop  the  average  stress  of 
1,000  psi  against  the  model.  In  the  initial  tests  of 
model  blocks  without  tunnels,  the  four  longitudinal 
hydraulic  jacks  were  driven  front  one  common  hand 
pump.  It  was  observed  that  any  tendency  of  the  longi¬ 
tudinal  loading  head  to  rotate  could  not  be  controlled 
with  this  system,  and  further  testing  was  done  with  an 
individual  hand  pump  for  each  jack,  ensuring  positive 
control  of  the  longitudinal  loading  head  movements. 


62 


Figure  Bl.V  So.  'malic  diagram  of  base  ami  loading  element  assembly. 


Figure  BI4.  Top  longitudinal  reaction  head. 
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